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ABSTRACT 


Recently  the  idea  of  deploying  non-collocated  sources  and  receivers  in  multistatic  sensor 
networks  (MSNs)  has  emerged  as  a  promising  area  of  opportunity  in  sonar  systems.  This 
thesis  addresses  point  coverage  sensing  problems  in  MSNs,  where  a  number  of  points  of 
interest  have  to  be  monitored  in  order  to  protect  them  from  hostile  underwater  assets.  We 
consider  discrete  “cookie  cutter”  sensors  as  well  as  various  diffuse  sensor  models.  By 
showing  that  the  convex  hull  spanned  by  the  targets  is  guaranteed  to  contain  optimal  sensor 
positions,  we  are  able  to  limit  the  solution  space.  Using  a  cookie  cutter  sensor  model, 
we  are  able  to  exclude  even  more  suboptimal  solutions  by  determining  range-of-the-day, 
source  and  receiver  circles. 

To  address  the  nonconvex  single-source  placement  problem,  we  develop  the  Divide 
Best  Sector  (DiBS)  algorithm,  which  quickly  provides  an  optimal  source  position  assuming 
fixed  receivers.  Starting  with  a  basic  implementation  of  DiBS,  we  show  how  incorporating 
advanced  sector  splitting  methods  and  termination  conditions  further  improve  the  algo¬ 
rithm.  We  also  discuss  two  ways  to  use  DiBS  to  find  multiple  source  positions  by  placing 
sensors  iteratively  or  simultaneously.  Finally,  we  conclude  that  DiBS  is  a  fast  and  simple 
algorithm  that  supports  a  wide  variety  of  sensor  models,  various  termination  conditions, 
and  objective  functions. 
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Executive  Summary 


Active  sonar  systems  have  long  constituted  an  important  sensing  mechanism  aboard  sub¬ 
marines  and  ships  in  anti-submarine  warfare  (ASW).  Recently,  however,  the  idea  of  deploy¬ 
ing  non-collocated  sources  and  receivers  has  emerged  as  a  promising  area  of  opportunity 
in  sonar  systems.  A  multistatic  sensor  network  (MSN)  consisting  of  a  number  of  non- 
collocated  sources  and  receivers  carries  a  number  of  advantages  such  as  more  complicated 
countermeasure  tactics,  improved  deployment  opportunities  as  well  as  reduced  costs  and 
noise  pollution.  The  complex  geometry  of  MSNs,  however,  makes  deployment  decisions 
like  number  of  sensors  and  their  locations  more  challenging. 

Out  of  three  sensing  problems  studied  in  the  literature  (area  search,  barrier  search,  and 
point  coverage)  this  thesis  addresses  point  coverage  problems  assuming  cookie  cutter  sen¬ 
sor  models  as  well  as  diffuse  sensor  models.  We  further  assume  stationary  targets  and  sen¬ 
sors  in  a  two-dimensional  Euclidean  plane  with  homogeneous  environment.  Additionally 
we  neglect  the  target’s  aspect  dependence  and  assume  independence  for  all  probabilities. 

We  start  by  showing  that  for  each  sensor  location  outside  the  convex  hull  encasing  all 
targets,  there  exists  a  position  inside  whose  detection  probability  is  at  least  as  good.  Based 
on  this,  we  limit  the  search  for  optimal  sensor  locations  to  the  convex  hull  of  the  targets. 

Next,  we  introduce  the  notion  of  range  of  the  day  circles  (RDCs)  for  cookie  cutter  sensor 
models,  i.e.,  circles  around  targets  whose  radii  are  the  range  of  the  day.  Utilizing  common 
points,  we  establish  the  minimal  set  of  clusters  of  RDCs,  G,  whose  cardinality  represents 
a  lower  bound  on  the  number  of  sensors  required  to  cover  all  targets.  Moreover,  G  lets  us 
define  an  upper  bound  on  the  number  of  detected  targets  with  a  restricted  number  of  sensors. 
Receiver  circles  (RCs)  and  source  circles  (SCs)  provide  additional  means  to  narrow  down 
possible  sensor  locations. 

Up  to  this  point,  we  can  formulate  integer  nonlinear  programs  (INLPs)  with  reasonable  run 
time  only  to  find  optimal  locations  for  exactly  one  source  and  one  receiver,  assuming  a 
cookie  cutter  sensor  model.  Hence  we  develop  the  Divide  Best  Sector  (DiBS)  algorithm  to 
find  a  single  optimal  source  position  assuming  fixed  target  and  receiver  positions  regardless 
of  the  applied  sensor  model.  The  algorithm  divides  the  area  of  possible  source  locations 


XV 


into  rectangular  sectors  and  evaluates  their  upper  bounds.  At  each  iteration  we  divide  the 
sector  with  the  highest  upper  bound  into  smaller  sector.  Eventually,  the  algorithm  reaches 
a  termination  condition  such  as  maximum  size  of  a  sector’s  longest  edge  or  maximum 
optimality  gap. 

We  proceed  to  investigate  and  assess  some  details  of  DiBS  starting  with  discovering  that  a 
small  number  of  newly  created  sectors  per  iteration  as  well  as  fewer  initial  sectors  greatly 
reduce  the  number  of  total  sectors  created  and  evaluated  by  the  algorithm.  This  subse¬ 
quently  leads  to  better  run  times. 

Another  method  to  further  improve  the  algorithm  is  rotating  the  plane  such  that  the  edges 
of  the  minimum-area  rectangle  encasing  all  targets  are  parallel  to  x-  and  y-  axis.  By  doing 
so,  we  minimize  the  sectors’  area  outside  the  convex  hull. 

Furthermore,  we  observe  that  uneven  rectangles  are  able  to  increase  the  total  number  of 
sectors  and  therefore  run  time.  We  discuss  two  ways  to  address  this  issue  by  adjusting  the 
termination  condition  and  by  creating  only  square  sectors. 

Finally,  we  explore  two  approaches  to  apply  DiBS  to  problems  that  require  finding  an 
optimal  placement  for  multiple  sources.  Here,  the  iterative  method  provides  a  lower  and 
upper  bound  assuming  a  cookie  cutter  sensor  model.  The  second  approach  places  sources 
simultaneously  and  evaluates  all  possible  combinations  of  sectors.  This  approach,  however, 
can  only  be  applied  to  small  numbers  of  sources. 
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CHAPTER  1: 
Introduction 


1.1  Background 

Active  sonar  systems  have  long  constituted  an  important  sensing  mechanism  aboard  sub¬ 
marines  and  ships  in  anti-submarine  warfare  (ASW).  In  a  typical  sonar  system,  a  ping  is 
sent  out  and  the  echo  yields  information  about  other  objects  in  the  area.  Recently,  however, 
the  idea  of  deploying  non-collocated  sources  and  receivers  has  emerged  as  a  promising  area 
of  opportunity  in  sonar  systems. 

A  multistatic  sensor  network  (MSN)  consisting  of  a  number  of  non-collocated  sources  and 
receivers  carries  a  number  of  advantages.  Cox  (1989,  p.  23)  states  that  “countermeasure 
tactics  are  greatly  complicated  if  the  target  does  not  know  the  position  of  the  receivers.” 
This  is  justified  by  the  fact  that  receivers  do  not  send  out  pings  and  thus  do  not  reveal 
their  locations.  Beyond  that,  procurement  estimates  of  the  United  States  Navy  (USN)  in¬ 
dicate  that  sources  might  cost  about  five  times  as  much  as  receivers  (USN,  2014).  Thus, 
deploying  more  receivers  than  sources  might  significantly  reduce  costs  without  sacrific¬ 
ing  performance.  Also,  a  “multistatic  system  can  employ  different  platforms  for  sources 
and  receivers.  A  ship  might  be  the  source,  while  the  receivers  are  sonobuoys”  (Washburn, 
2010,  p.  1).  On  top  of  this,  a  ping  might  be  received  by  multiple  receivers.  Coon  (1997) 
and  Simakov  (2008)  discuss  how  to  merge  multiple  detections  into  a  single  alert  that  is 
more  precise  and  eliminates  some  of  the  false  alarms  that  occur  on  traditional  (monostatic) 
sonar  systems. 

The  performance  of  a  bistatic  sonar  system,  i.e.,  an  MSN  with  exactly  one  source  and  one 
receiver,  is  significantly  more  difficult  to  model  than  a  monostatic  sonar  system,  leading 
to  challenges  in  optimal  deployment  and  usage.  The  primary  source  of  these  challenges 
lies  in  the  differences  of  the  geometry  of  both  systems.  In  a  monostatic  sonar  system  the 
detection  probability  is  mainly  related  to  the  distance  between  the  sonar  device  and  the  po¬ 
tential  target.  This  relationship  is  more  complicated  in  a  bistatic  model,  where  the  detection 
probability  depends  on  the  product  of  the  target-source  and  target-receiver  distances  (Cox, 
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Figure  1.1:  Geometry  of  a  Bistatic  Sonar  System  -  The  detection  probability  for  target  t 
depends  on  the  product  of  the  target-source  and  target-receiver  distance. 

1989)  as  displayed  in  Figure  1.1.  Moreover,  the  analytical  challenges  are  exacerbated  by 
MSNs  involving  multiple  sources  and  receivers. 

1.2  Literature  Review 

Researching  sonar  technology  and  MSNs  in  particular  is  a  wide  field  with  most  diverse 
subcategories.  We  will  discuss  a  selection  in  this  section,  starting  with  the  geometry  of 
multistatic  models  and  following  with  examples  for  usage  and  deployment. 

1.2.1  The  Geometry  of  Multistatic  Sensor  Networks 

Multiple  authors  describe  the  geometry  of  MSNs.  Cox  (1989)  analyzes  the  relationship  be¬ 
tween  monostatic  and  bistatic  active  sonars.  He  derives  that  a  detection  probability  contour, 
i.e.,  a  contour  consisting  of  all  locations  for  a  target  t  with  the  same  detection  probability, 
is  defined  by  the  constant  product 


(1.1) 


where  pt,s,r  is  constant,  and  dt^s  and  dt^r  are  the  distances  from  a  target  t  to  a  source  5 
and  to  a  receiver  r  respectively.  The  distances  are  also  illustrated  in  Figure  1.1.  Those 
contours  are  geometric  figures  known  as  Cassini  ovals  shown  in  Figure  1.2.  The  team  of 
Matthew  Fewell  and  Sylvia  Ozols  from  the  Australian  Defence  Science  and  Technology 
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^source;  A  receiver 


Figure  1.2:  Examples  of  Cassini  Ovals  -  Cassini  ovals  are  shown  with  respect  to  source  and 
receiver  distance  ds.r  where  po  denotes  the  range  of  the  day.  If  source  and  receiver  are  collocated 
then  the  Cassini  oval  is  a  circle  that  corresponds  to  a  monostatic  sensor.  Increasing  the  distance 
changes  the  shape  from  oval  over  dog-bone  shape  up  to  two  separate  egg  shapes. 


Organisation  (DSTO)  published  multiple  reports  on  MSNs.  Fewell  and  Ozols  (2011,  p.  4) 
state  that  the  constant  pt,s,r  from  Equation  (1.1)  is  the  equivalent  range  to  the  range  from  a 
monostatic  sonar  system,  which  enables  us  to  compute  detection  probabilities  for  bistatic 
sensor  networks.  The  authors  point  out  some  issues  with  this  model,  e.g.,  “that  not  only  [the 
detection  probability  is]  high  at  a  receiver  but  also  the  value  is  unaffected  by  the  source- 
receiver  distance.  This  seems  inadequate  (Fewell  &  Ozols,  2011,  p.  7).”  They  list  various 
possibilities  to  modify  the  model  to  account  for  these  issues. 

There  are  multiple  ways  to  determine  Pt,s,r,  the  probability  to  detect  target  t  with  source  5 
and  receiver  r  in  a  multistatic  model  as  displayed  in  Figure  1.3.  The  simplest  one  is  the 
cookie  cutter  sensor  model,  also  known  as  the  definite  range  sensor  in  literature.  Based  on 
the  range  of  the  day  po,  i.e.,  the  distance  from  a  monostatic  sensor  to  a  target  where  the 
detection  probability  is  50%,  we  define  Pt^s,r  as 

_  Jl  if  pb  <  Pt,s,r  <  Po, 

[0  otherwise. 

Other  names  for  the  range  of  the  day  are  the  detection  range  or  i?5o.  The  blind  zone  pt, 
is  the  area  where  targets  cannot  be  detected  because  their  echoes  arrive  at  nearly  the  same 
time  as  the  ping  from  the  source.  Then  the  total  detection  probability  for  target  t,  Pt,  i.e.. 
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Figure  1.3:  Sensor  Models  -  The  probability  curves  for  the  three  sensor  models  cookie  cutter, 
Fermi  and  exponential  function  are  displayed.  Range  is  expressed  as  multiples  of  the  range  of 
the  day  po.  All  models  have  a  Pt^s,r  =  0.5  at  Pt,s,r/po  =  1  to  be  consistent  with  the  definition  of 
the  range  of  the  day. 


the  probability  to  detect  target  t  with  all  sensors  in  scene,  is  defined  as 


Pt=  max  Pt^s,r-  (1-3) 

{s,r)eSxR 


This  means,  if  at  least  one  pair  of  sensors  detects  target  t  then  Pt  =  I,  otherwise  Pt  =  0. 
Although  Equations  (1.2)  and  (1.3)  form  an  analytically  convenient  model,  it  lacks  some 
features  of  real  sensors  such  as  a  gradually  decreasing  Pt^s,r  with  increasing  p. 


Hence,  the  DSTO  team  proposes  two  diffuse  sensor  models  in  (Fewell  &  Ozols,  2011). 
The  first  is  the  Fermi  function^ 


Ptj:.r  — 


{ _ ^ _ 

J  1  -I-  \Q{.Pt,s,r/P0~'^)/b 

1° 


if  Pt,s,r  ^  Pb^ 

Otherwise. 


(1.4) 


The  diffusivity  parameter  b  determines  how  rapidly  probability  values  change  when  chang¬ 
ing  range  p.  As  b  0,  the  Fermi  function  approaches  the  cookie  cutter  model. 


'Discovered  in  1926  by  Enrico  Fermi  and  Paul  Dirac  when  researching  electron  behavior. 


4 


The  second  model  is  the  exponential  function 


10^0.30103p(„5^r/ Po 

0 


if  Pt,s,r  ^  Pb: 

Otherwise. 


(1.5) 


The  exponent  is  chosen  such  that  if  pt,s,r/ Po  =  1  the  detection  probability  =  0.5.  This 
is  important  to  be  consistent  with  the  definition  of  the  range  of  the  day. 

For  diffuse  sensor  models  we  need  a  different  approach  to  evaluate  the  total  detection  prob¬ 
ability.  Thus,  we  define  a  more  general  version  of  Equation  (1.3)  to  compute  Pt  regardless 
of  the  chosen  sensor  model  assuming  all  probabilities  are  independent  as 


Pt=^-  n  (1-6) 

{s,r)eSxR 

In  order  to  account  for  different  weighted  targets,  we  introduce  Vt  as  the  value  of  target 
t.  Hence  we  denote  the  expected  reward  for  target  t  as  VtPt-  If  all  targets  have  the  same 
value,  we  simplify  the  expected  reward  for  target  t  by  using  Pf.  Based  on  this,  we  can 
define  multiple  objective  functions  depending  on  the  chosen  sensor  model.  While  maxi¬ 
mizing  the  total  expected  reward  is  our  only  objective  when  using  a  cookie  cutter  model, 
we  can  also  choose  to  maximize  the  minimum  expected  reward.  The  objective  functions 
are  summarized  in  Table  1.1. 


Objective 

Sensor  Model 

Formula 

maximize  total  expected  reward 

cookie  cutter 

max  Y,  VtPt 
teT 

maximize  average  expected  reward 

diffuse 

max  Y  VtPt/\T\ 

t^T 

maximize  minimum  expected  reward 

diffuse 

max  min  V;/) 

teT 

Table  1.1:  Objective  functions  -  Three  different  objective  functions  are  possible  depending  on 
the  chosen  sensor  model,  where  V;  denotes  the  value  or  weight  of  target  t. 


1.2.2  Deployment  and  Usage 

There  are  different  approaches  to  quantify  the  effectiveness  of  the  deployment  and  usage 
of  multistatic  sensors.  A  strategic  rather  than  a  tactical  approach  is  analyzed  by  Washburn 
(2010).  He  assumes  that  sources  and  receivers  are  deployed  uniformly  at  random  within 
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some  region.  The  author  develops  methods  for  approximating  the  detection  probabilities 
as  a  function  of  number  of  sensors.  The  most  interesting  aspect  about  this  approach  is  the 
fact  that  it  does  not  need  to  consider  the  geometric  arrangement  of  the  sensors. 

The  DSTO  analyzes  multiple  scenarios.  Fewell  and  Ozols  (2011)  compare  a  field  of  mono¬ 
static  sonar  systems  with  that  of  a  field  of  similar  sonars  operated  multistatically,  where 
sources  and  receivers  are  collocated.  The  direct  comparison  of  the  two  modes  of  operation 
reveals  that  the  correct  choice  of  sensor  models  affects  the  outcome.  Using  a  cookie  cutter 
sensor  model,  the  researchers  find  no  advantage  to  the  MSN.  However,  using  the  expo¬ 
nential  model,  the  researchers  achieved  the  same  results  in  an  MSN  with  about  one  quarter 
the  number  of  sensors  as  the  monostatic  sensors.  This  finding  is  important  since  it  reduces 
the  number  of  pings  in  a  given  field.  As  soon  as  a  ping  is  sent  out,  a  hostile  submarine 
knows  the  location  of  the  source  and  will  depart,  which  consequently  makes  it  harder  to 
detect.  Another  reason  to  reduce  the  number  of  pings  is  the  artificial  stress  for  sea  dwellers 
produced  by  sonar  systems. 

There  are  three  main  types  of  sensing  coverage  problems  in  the  literature:  area  search, 
barrier  search,  and  point  coverage.  Using  the  same  sensor  models  as  in  (Fewell  &  Ozols, 
2011),  Ozols  and  Fewell  (2011)  compare  various  layouts  for  sensor  positions  in  a  MSN, 
where  a  large  area  has  to  be  covered.  From  the  27  layouts  tested,  four  are  recommended  for 
use  depending  on  the  ratio  of  source  and  receiver  cost.  Additionally,  in  a  classified  report, 
Fewell,  Ozols,  and  Rzetelski  (2011)  discuss  the  deployment  of  MSNs  for  a  barrier  search 
problem.  Gong,  Zhang,  Cochran,  and  Xing  (2013)  address  the  barrier  search  problem  by 
finding  the  optimal  deployment  of  multistatic  sensors  to  maximize  the  worst-case  intrusion 
detectability.  The  authors  show  that  a  balanced  structure  yields  the  optimal  solution. 

Craparo  and  Karata§  (2014)  engage  the  point  coverage  problem  for  MSNs.  The  idea  behind 
is  to  monitor  a  number  of  points  of  interest  (Pols),  e.g.,  oil  platforms,  aircraft  carriers, 
ports,  etc.,  in  order  to  detect  hostile  underwater  assets.  The  authors  assume  fixed  Pols 
and  receivers,  and  discuss  various  approaches  to  optimally  place  multiple  sources.  To  our 
knowledge  this  is  the  only  existing  study  considering  the  point  coverage  problem  for  MSNs. 
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1.3  Objectives 

From  the  three  main  types  of  sensing  coverage  problems  (barrier  search,  area  search  and 
point  coverage),  this  thesis  focuses  on  the  point  coverage  problem  for  MSNs.  The  overall 
question  is:  where  do  we  place  sources  and  receivers  such  that  the  Pols  are  optimally 
covered?  To  be  in  accordance  with  the  terminology  of  other  literature,  ‘Pols’  are  called 
‘targets’,  and  instead  of  ‘covering’  we  say  ‘detecting’.  To  engage  this  problem,  we  consider 
the  following  questions. 

•  What  are  basic  observations  regarding  point  coverage  sensing  in  MSNs?  Particularly, 

-  Is  the  convex  hull  spanned  by  the  targets  guaranteed  to  contain  the  optimal 
locations  of  sources  and  receivers? 

-  Are  there  other  ways  to  exclude  suboptimal  solutions? 

•  How  can  we  find  an  optimal  placement  in  MSNs?  Particularly, 

-  Are  there  models  or  algorithms  that  represent  point  coverage  sensing? 

-  How  do  the  developed  means  perform  and  how  can  they  be  improved? 


1.4  Scope,  Limitations,  and  Assumptions 

For  the  developed  models  we  assume  that  all  targets  and  sensors  are  stationary  and  exist 
in  a  two-dimensional  Euclidean  plane  with  homogeneous  environmental  conditions.  In  the 
course  of  this  study,  all  three  detection  probabilities  shown  in  Figure  1.3  are  considered. 
When  we  build  our  models,  however,  we  start  with  the  simple  cookie  cutter  model  and  later 
advance  to  the  diffuse  sensor  models.  The  effect  of  the  blind  zone  can  be  greatly  reduced 
by  pulse  compression  as  shown  by  Fewell  and  Ozols  (2011,  p.  11).  For  this  reason,  we  do 
not  consider  blind  zones  in  this  thesis,  i.e.,  we  assume  pb  =  0.  Furthermore,  we  neglect  the 
target’s  aspect  dependence  and  assume  all  probabilities  are  independent.  A  set  of  notations 
used  throughout  the  study  is  summarized  in  Table  1.2. 

1.5  Contributions  And  Outline 

In  the  first  part  of  this  thesis  we  explore  some  general  observations  about  point  coverage 
sensing  with  MSNs.  We  prove  that  a  convex  hull  encasing  the  targets  contains  optimal 
sensor  locations.  Furthermore,  we  show  the  importance  of  various  circles  and  the  clusters 
of  targets  induced  by  their  intersections  for  cookie  cutter  sensor  models.  We  show  that  these 
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Element 

Set 

Deseription 

t 

T 

target 

s 

S 

souree 

r 

R 

reeeiver 

Vt 

M+ 

value  of  target  t 

{xi,yi) 

eoordinates  of  objeet  /,  V/  G  T  U  5  U/? 

M+ 

Euelidean  distanee  between  objeet  i  and  j,  \/i,j  ^  TUSUR 

Po 

M+ 

range  of  the  day 

Pt,s,r 

M+ 

equivalent  monostatie  range  for  target  t,  souree  s,  and  reeeiver  r 

Pt,s,r 

M+ 

probability  to  deteet  target  t  with  souree  s  and  reeeiver  r;  Pt^s,r  <  1 

P, 

M+ 

total  deteetion  probability  for  target  t;  P?  <  1 

Table  1.2:  Notation  -  This  table  summarizes  the  notation  used  throughout  the  study. 


clusters  can  be  used  to  define  bounds  on  the  sensing  problem.  Additionally,  we  discuss  how 
integer  nonlinear  programs  (INLPs)  can  be  used  to  find  an  optimal  solution. 

The  second  part  focuses  on  the  development  and  enhancement  of  a  new  algorithm  to  find 
the  optimal  source  position,  assuming  fixed  targets  and  receivers  in  place.  We  investigate 
many  details  of  the  algorithm  and  assess  methods  and  means  to  improve  it. 
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CHAPTER  2: 

Observations  On  Point  Coverage  Sensing 


This  chapter  analyzes  the  point  coverage  scenario  and  discovers  some  basic  observations. 
Those  observations  help  to  limit  the  solution  space  and  exclude  infeasible  or  non-optimal 
solutions.  We  also  show  how  to  apply  the  gained  insights  and  whether  the  outcome  justifies 
utilizing  the  developed  method. 

2.1  The  Convex  Hull 

A  simple  approach  to  limit  the  set  of  potential  sensor  locations  would  be  defining  a  rectan¬ 
gle  that  contains  all  Pols.  Even  though  this  can  easily  be  implemented,  we  want  to  exclude 
as  many  locations  as  possible.  A  more  sophisticated  approach  involves  looking  at  the  con¬ 
vex  hull  spanned  by  the  targets.  As  vividly  described  in  (De  Berg,  Van  Kreveld,  Overmars, 
&  Schwarzkopf,  2000,  p.  3),  one  can  imagine  the  targets  as  nails  sticking  out  of  the  plane. 
If  we  hold  an  elastic  rubber  band  around  the  nails  and  let  it  go,  the  enclosed  area  will  be 
the  convex  hull. 

2.1.1  Properties  Of  The  Convex  Hull 

The  example  in  Figure  2.1  shows  a  convex  hull  for  a  particular  set  of  targets.  It  also  shows 
that  all  targets  are  either  inside  the  hull  or  at  one  corner  of  the  resulting  polygon.  The  latter 
are  called  vertices.  Hence  we  define  the  set  C  as 

C  =  {tGr|tisa  vertex  of  the  convex  hull  of  T}. 

The  convex  hull  is  defined  as  the  set  Conv(r),  that  contains  all  points,  that  are  spanned 
by  the  vertices  in  T .  It  is  easy  to  see  that  Conv(r)  =  Conv(C).  If  all  nails  that  are  not 
vertices  are  removed  in  the  previously  mentioned  example,  the  rubber  band  would  still 
enclose  the  same  area.  Alongside  this  picturesque  description  of  a  convex  hull,  there  also 
exists  a  mathematical  definition.  Each  point  p  G  Conv(C)  with  coordinates  {xp,yp)  can  be 
written  as  a  linear  combination  of  the  vertices  in  C.  This  is  shown  in  Equations  (2.1),  where 
0  <  Ar  <  1  for  all  t  G  C. 


9 


78 


Figure  2.1:  Convex  Hull  Example  -  A  Convex  Hull  spanned  by  a  set  of  targets.  Its  vertices  are 
Ti,  T(i,  Tg,  Tio  and  Ti. 


X 


P  ~ 


teC 


yp  = 

leC 


1  =  IA,. 

teC 


(2.1) 


Theorem  2.1  deseribes  the  relationship  between  the  eonvex  hull  and  the  deteetion  proba¬ 
bility  for  sensor  p,  where  Pt{p)  denotes  if  p  is  a  souree  and  Pt^s,p  otherwise. 

Theorem  2.1.  For  every  sensor  position  p  ^  Conv(C),  there  exists  a  position  p'  G  Conv(C), 
such  that  Pt{p)  <  Pt{p')yt  G  T. 


Proof.  Assume  sensor  position  p  ^  Conv(C).  Let  p'  G  Conv(C)  denote  the  position  inside 
the  eonvex  hull  with  the  shortest  distanee  to  p.  Then  p'  is  either  on  an  edge  of  the  eonvex 
hull  or  eolloeated  with  one  of  its  vertiees  as  shown  in  Figure  2.2.  In  both  oases,  the  eonvex 
hull  and  with  it  all  targets  t  G  T  are  behind  an  imaginary  line  perpendioular  to  the  line  from 
p  to  p'  orossing  that  line  at  p' .  Otherwise  there  would  be  another  p'  that  is  eloser  to  p, 
whioh  is  a  oontradiotion. 
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Figure  2.2:  Theorem  2.1  Proof  -  There  exist  two  cases  for  the  closest  position  p'  e  Conv(C)  to 
p.  The  left  figure  displays  p'  on  an  edge  of  the  convex  hull,  while  p'  is  collocated  with  a  vertex 
on  the  right.  In  both  cases  all  targets  f  e  T  are  right  of  the  dashed  line,  which  is  perpendicular 
to  the  line  from  p  to  p'  and  crosses  p' . 


Let  t  be  a  target  inside  the  convex  hull.  Without  loss  of  generality,  we  assume  that  p  and 
p'dSQ  lying  on  a  line  parallel  to  the  x-axis.  Then  the  horizontal  distance  between  p  and  t  is 
\xp  —Xt\  =  \xp—  Xpf I  +  \xp/  —  Xt I  and  we  derive 

dip  =  {^p-^tf  +  {yp-ytf 

=  {\xp-Xpi\  +  \xpi -Xt\f +  {yp-ytf 

>  [Xp!  -Xtf  +  {ypi  -Jtf 


Let  pt{p)  denote  p?,p,r  and  pt^s,p  if  p  is  a  source  or  receiver  respectively.  Thus  the  equivalent 
monostatic  range  for  p'  is 

Pt{p')  =  \Jdt^p’dt^x  <  \/dt,pdt^x  =  Ptip), 

where  x  is  the  second  type  of  sensor  needed  in  an  MSN.  It  follows  that  the  detection 
probability  for  a  target  t  using  sensor  p' ,  Pt{p')  <  Pt{p)  since  all  sensor  model  functions 
are  monotonically  nonincreasing  with  distance  as  is  shown  in  Figure  1.3.  □ 


We  conclude  that  applying  Theorem  2.1  we  are  able  to  place  all  sources  and  receivers 
inside  the  convex  hull  without  sacrificing  detection  probabilities.  Hence,  we  are  now  able 
to  constrain  our  solution  space,  i.e.,  we  can  find  optimal  positions  for  each  sensor  inside 
the  convex  hull.  We  now  describe  and  assess  two  ways  to  construct  the  convex  hull. 
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2.1.2  A  Linear  Program  To  Find  Vertices 


The  first  method  uses  a  linear  program  (LP)  to  find  the  vertices  of  the  convex  hull.  We  use 
the  fact  that  vertices  themselves  cannot  be  written  as  a  linear  combination  of  other  points, 
whereas  targets  that  are  not  vertices  can.  Based  on  this  insight,  we  can  define  the  following 
LP. 

Indices  and  Sets: 

t,t’ &T  targets. 


Data  [units]: 

Xt  X  coordinate  of  target  t  [unitless], 

yt  y  coordinate  of  target  t  [unitless]. 

Decision  Variables: 

A/  coefficient  for  target  t's  position  subject  to  target  t. 

Formulation: 


min 

A 


teT 


Subject  to 

Xt'  =  ^ 
teT 

Vt'  e  T, 

(2.2a) 

yf  =  E  ^t'yt 

teT 

'it'  e  T, 

(2.2b) 

teT 

it'  e  r, 

(2.2c) 

0  <  A/' 

itj'  G  T. 

The  LP  tries  to  minimize  the  number  of  target  locations  not  expressed  as  a  convex  combi¬ 
nation  of  other  target  locations.  For  a  non-vertex  target  t  A/  can  be  zero,  since  t’s  position 
can  be  represented  by  a  linear  combination  of  the  vertices.  The  only  possible  value  for 
vertices,  however,  is  one.  Thus,  we  can  define  the  set  of  vertices  C  as 

C  =  {teT\Xl  =  \}. 
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Implementing  and  solving  this  LP  in  General  Algebraic  Modeling  System  (GAMS)  runs 
very  quickly  for  small  numbers  of  targets.  Figure  2.3,  however,  shows  that  more  targets 
result  in  an  disproportional  increase  in  run  times  when  using  the  simplex  algorithm.  In 
fact,  Klee  and  Minty  (1972)  show  that  in  a  worst-case  scenario,  the  simplex  algorithm  can 
reach  exponential  complexity. 


Figure  2.3:  Convex  Hull  LP  Run  Time  -  For  small  numbers  of  targets  the  LP  runs  very  quickly. 
Run  Times,  however,  grow  very  quickly  with  the  number  of  targets. 


2.1.3  The  Graham  Scan  Algorithm 

Another  algorithm  to  find  the  vertices  of  a  convex  hull  with  complexity  0(|r|log|r|)  is 
the  Graham  scan  described  in  (Graham,  1972).  The  algorithm  starts  with  finding  tQ  G  T , 
the  target  with  the  smallest  value  for  yt.  If  there  are  multiple  targets  that  share  the  smallest 
yt,  we  pick  the  one  with  the  smallest  Xt  out  of  the  candidates. 

Next  we  compute  the  angles  Of  each  target  t  G  T  makes  with  tQ  and  the  x-axis  using  the 
formula 


0t  =  atan^ {yt  - yt^ .Xt-Xt^).  (2.3) 

Sorting  the  targets  by  9t  starting  with  the  smallest,  the  algorithm  now  considers  each  target 
as  possible  vertex.  At  each  step  it  is  determined  whether  the  current  target  t  and  its  two 
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predecessors  t  —I  and  t  —  2  make  a  turn  clockwise  or  counter  clockwise  by  calculating 


c  =  {xt^i  -xt-2)  X  (yt-yt-i)  -  (yt-i  -yt-2)  x  {xt-xt-2).  (2.4) 

If  the  result  is  positive  the  points  make  a  counter  clockwise  turn  and  the  next  target  in  the 
ordered  list  is  considered.  If  c  <  0  then  the  turn  is  counter  clockwise  and  the  middle  target 
t  —  1  is  identified  as  non- vertex.  At  c  =  0  all  three  targets  are  collinear,  which  also  leads  to 
discarding  t  Figure  2.4  displays  the  working  of  the  Graham  scan  algorithm. 


Figure  2.4:  Graham  Scan  -  The  working  of  the  Graham  scan  is  displayed.  Targets  are  considered 
in  order  of  the  angle  they  make  between  the  lowest  target  and  the  x-axis.  In  this  example  at  the 
second  step  the  algorithm  detects  a  clockwise  turn,  which  leads  to  discarding  the  middle  target 
on  step  3. 


The  low  complexity  of  the  Graham  scan  is  a  result  of  the  simple  formula  in  Equation  (2.4) 
to  decide  the  complex  issue  about  clockwise  or  counterclockwise  turns.  Outputting  the 
vertices  of  the  convex  hull  in  a  counter  clockwise  order  is  another  advantage  we  later  exploit 
in  Section  3.4.3. 


2.2  Range  Of  The  Day  Circles 

This  section  analyzes  the  relationship  between  the  range  of  the  day,  po  and  the  cookie  cutter 
sensor  model.  We  already  know  from  Equation  (1.2)  that  a  target  is  detected  if  and  only 
if  its  equivalent  monostatic  range  pt,s,r  <  Po-  Besides  that,  we  introduce  range  of  the  day 
circles  (RDCs),  i.e.,  circles  around  targets  with  radius  po  as  shown  in  Eigure  2.5.  In  the 
following,  we  discuss  methods  that  use  RDCs  that  bound  the  problem.  Even  though  we 
develop  these  methods  for  cookie  cutter  sensor  models,  they  also  apply  to  problems  where 
we  have  to  meet  a  particular  detection  probability,  e.g.,  80%.  Here  we  simply  change  the 
radius  to  a  value  that  corresponds  with  the  demanded  probaility. 
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Figure  2.5:  Range  Of  The  Day  Circles  Example  -  Assuming  a  cookie  cutter  sensor  model, 
range  of  the  day  circles  determine  bounds  for  the  problem.  Each  circle  has  a  target  for  center 
and  po  for  radius. 


2.2.1  Properties 

An  important  observation  is  the  relationship  between  sensor  positions  relative  to  RDCs  and 
a  target’s  detection  probability. 

Theorem  2.2.  A  target  t  is  detected  only  if  at  least  one  sensor  is  inside  t  ’s  RDC. 


Proof.  To  arrive  at  a  contradiction,  we  assume  all  sensors  are  located  outside  the  RDC  of 
target  t  and  Pt^s,r  =  1  for  some  5  G  5  and  r  eR.  Then  dt^s  >  Po,  V5  G  S  and  dt^r  >  Po,  Vr  G  R. 
It  follows  that 


Pt,s,r  —  y/ dt^s  X  dt^r 

>  a/Po  X  po  =  po- 

But  then  according  to  Equation  (1.2)  Pt^s,r  =  0,  which  is  a  contradiction.  □ 

This,  however,  does  not  mean  that  we  can  find  an  optimal  sensor  placement  by  solely  look- 
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Figure  2.6:  Range  Of  The  Day  Circles  Properties  -  The  left  plot  illustrates  that  an  optimal 
placement  can  be  found  by  placing  the  source  outside  RDCs.  Forcing  the  inside  the  circles 
results  in  losing  coverage  of  targets.  The  right  plot  shows  an  example  where  an  optimal  receiver 
position  is  outside  RDC  intersections.  Here,  it  is  not  possible  to  cover  all  targets  by  placing  all 
receivers  into  RDC  intersections. 


ing  at  positions  inside  RDCs.  The  left  plot  in  Figure  2.6  demonstrates  a  eounterexample 
where  the  optimal  position  for  the  souree  is  loeated  between  both  eliques  of  targets  out¬ 
side  any  RDCs.  Moving  this  souree  inside  a  RDC  of  one  of  the  eliques  results  in  losing 
eoverage  of  targets  in  the  other  elique. 

A  reasonable  heuristie  may  involve  ehoosing  to  plaee  reeeivers  into  the  eireles  and  let  the 
sourees  be  the  eonneeting  elements  of  the  MSN  in  order  to  keep  the  number  of  sourees 
as  small  as  possible.  However,  the  impression  that  optimal  reeeiver  positions  ean  only 
be  found  in  RDC  interseetions,  unless  a  eirele  does  not  have  an  interseetion  with  another 
eirele,  is  fallaeious  as  pietured  by  the  right  plot  in  Figure  2.6.  Reeeiver  Ri  eannotbe  moved 
into  the  interseetion  of  Ti’s  and  T2S  RDCs  without  saerifioing  eoverage  of  T\.  Covering 
targets  Tt,  and  inhibits  the  souree  from  moving  eloser  to  T\ .  We  take  a  eloser  look  at  this 
example  in  Seetion  2.2.3. 


2.2.2  Clusters 

Nonetheless  we  ean  use  RDCs  and  their  interseetions  to  form  elusters.  A  eluster  G  C  T 
is  a  maximal  set  of  targets,  where  the  RDCs  of  all  targets  t  G  G  have  at  least  one  point  in 
eommon.  Sinee  we  require  a  eluster  to  be  maximal,  it  is  not  possible  to  add  another  target 
to  a  eluster.  In  our  example  from  Figure  2.5  we  ean  find  the  elusters 


G2  =  {r3,r7},  G3  =  {r2,r6,r9}, 

G5  =  {T4,Ts}- 


Gi  =  {Ti}, 

G4  =  {r4,r5,rio}, 
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We  make  some  observations  about  clusters.  First  of  all,  each  target  t  G  T  occurs  in  at  least 
one  cluster.  If  there  are  no  intersections  with  circles  from  other  targets,  a  target  forms  its 
own  cluster,  e.g.,  Ti  is  the  only  element  of  G\.  Additionally,  it  is  possible  that  a  target  is  a 
member  of  multiple  clusters,  e.g.,  T/^  G  G4  and  Moreover,  it  is  important  to  notice 

that  even  though  at  a  first  glance  it  looks  like  the  clique  of  targets  on  top  of  Figure  2.5  (r4, 
T^,  Tg  and  Tiq)  all  belong  together,  they  actually  form  two  separate  clusters  since  Tg  is  only 
connected  to  r4. 

Figure  2.7  shows  a  special  case:  the  left  plot  forms  exactly  one  cluster,  while  the  right  plot 
forms  three  clusters  {ri,r2},  {r2,r3}  and  {ri,r3}  since  the  three  targets  do  not  have  a 
single  point  in  common,  even  though  the  targets’  RDCs  mutually  intersect. 


Figure  2.7:  Clusters  Special  Case  -  The  targets  in  the  left  plot  form  exactly  one  cluster.  The 
targets  on  the  right,  however,  form  the  three  clusters  {ri,r2},  {T2^T^}  and  since  the 

three  targets  do  not  have  a  single  point  in  common. 

Based  on  this,  we  define  the  minimal  set  of  clusters  G  as  the  smallest  set  of  clusters  that 
contains  all  targets  t  G  T.  For  the  right  plot  in  Figure  2.7  a  minimal  set  of  clusters  is 
G  =  {{Ti,  r2},  {r2,  r3}}.  Furthermore,  we  formulate  the  following  theorem. 

Theorem  2.3.  A  lower  bound  on  the  number  of  sensors  required  to  detect  all  targets  is  |G|. 

Proof.  To  arrive  at  a  contradiction,  we  assume  that  all  targets  are  covered  with  |G|  —  1 
sensors.  It  follows  from  Theorem  2.2  that  at  least  one  sensor  has  to  be  inside  each  target’s 
RDC.  Since  the  circles  of  each  cluster  have  at  least  one  point  in  common,  putting  a  sensor 
on  that  point  results  in  having  a  sensor  in  each  RDC  of  this  cluster.  Moreover,  since  clusters 
are  maximal,  we  cannot  add  another  target.  Hence,  one  cluster  and  therefore  at  least  one 
target  remains  without  a  sensor  in  its  respective  RDC.  □ 
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As  a  result  from  Theorem  2.3  we  ean  define  |G|  as  a  lower  bound  for  the  number  of  sensors 
required  to  deteet  all  targets.  By  the  same  token,  we  are  able  to  find  an  upper  bound  for  the 
number  of  targets  we  ean  eover  with  a  given  number  of  sensors  n  by  ehoosing  n  sets  from 
G,  sueh  that  the  number  of  eontained  targets  is  maximized. 

2.2.3  Other  Circles 

Next  to  RDCs,  there  exist  other  eireles  that  provide  insight  into  where  to  plaee  either 
sourees  or  reeeivers.  Craparo  and  Karata§  (2014)  introduce  the  notion  of  the  detection 
disk  for  each  target  t  and  receiver  r  that  has  t  for  center  and  its  radius  defined  by 


In  order  to  better  distinguish  detection  disks  from  a  second  type  of  circles  we  introduce  in 
this  section,  we  refer  to  them  as  receiver  circles  (RCs).  In  order  to  detect  a  target  t,  we  have 
to  place  a  source  inside  fs  RC. 

Figure  2.8  shows  the  example  from  Figure  2.6  augmented  with  its  RCs.  Here  the  red  and 
blue  circles  are  RCs  for  Ri  and  R2,  respectively.  Using  Ri  and  one  source  we  are  able  to 
detect  Ti  and  T2  simultaneously  since  their  RCs  overlap.  If  we  want  to  detect  Tt,  or  T4, 
however,  we  are  not  able  to  detect  more  than  one  target  since  their  Ri  RCs  do  not  overlap 
with  other  Ri  RCs.  On  the  other  hand,  the  R2  RCs  of  and  T4  overlap  and  even  have  a 
common  point  with  the  RC  of  T2,  such  that  placing  a  source  here  results  in  detecting  T2,  Tt, 
and  T4  with  Ri. 

By  the  same  token,  we  can  define  source  circles  (SCs)  whose  radius  is  determined  by  a 
source  location  as 


Figure  2.9  displays  the  SCs  in  our  previous  example.  With  the  shown  fixed  source  position, 
placing  a  receiver  at  the  intersection  of  T^’s  and  r4’s  SCs  ensures  detecting  both  targets  as 
well  as  T2.  Placing  a  second  receiver  somewhere  into  Ti’s  SC  always  detects  T2  also  since 
T\s  circle  is  completely  overlapped  by  r2’s  SC. 
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Figure  2.8:  Receiver  Circles  Example  -  Red  and  blue  circles  represent  receivers  i?i's  and  7?2's 
receiver  circles  respectively.  Placing  a  source  at  RC  intersections  ensures  detection  of  respective 
targets. 


Figure  2.9:  Source  Circles  Example  -  Placing  a  receiver  at  the  intersection  of  targets  T^'s  and 
T4’s  source  circles  guarantees  detecting  both  targets.  A  receiver  somewhere  in  Ti's  SC  always 
detects  T2  as  well  because  of  the  total  overlapping. 


Assuming  the  positions  of  one  type  of  sensors  is  already  fixed,  utilizing  RCs  and  SCs  is  a 
way  to  decide  positions  for  the  remaining  type.  Nevertheless,  a  large  number  of  targets  and 
sensors  makes  their  manual  use  hardly  practical. 
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2.3  Utilizing  Nonlinear  Programs 

We  now  investigate  the  usage  of  INLPs  to  find  optimal  sensor  placements  in  MSNs.  The 
first  model  finds  an  optimal  solution  for  a  bistatic  sensor  network,  i.e.,  exactly  one  source 
and  one  receiver,  assuming  a  cookie  cutter  sensor  model.  Then  we  discuss  how  INLPs 
perform  on  extended  networks. 


2.3.1  Bistatic  Cookie  Cutter  Sensor  Networks 

For  this  model  we  start  from  the  premise  that  we  already  found  the  set  of  vertices  C,  using 
the  Graham  scan,  for  example.  Then  we  formulate  the  following  INLP  to  find  optimal 
positions  for  a  bistatic  sensor  network  assuming  a  cookie  cutter  sensor  model. 

Indices  and  Sets: 

t  G  r  targets, 

p  &  P  =  {s,r}  type  of  sensor,  5  for  source  and  r  for  receiver. 

Data  [units]: 

Xt  X  coordinate  of  target  t  [unitless] , 

yt  y  coordinate  of  target  t  [unitless], 

Vt  weight  of  target  t  [unitless], 
po  range  of  the  day  [miles], 

C  vertex  set,  C  =  {t  G  r|t  is  a  vertex}  [unitless], 

Mt  penalty  for  not  covering  target  t  [unitless]. 

Decision  Variables: 

Xj’  coefficient  for  sensor  p’s  position  subject  to  vertex  t, 

Xp  X  coordinate  of  sensor  p, 
yp  y  coordinate  of  sensor  p, 
ht  binary,  1  if  target  t  is  not  detected,  0  otherwise. 
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Formulation: 


min  Y^Vfht 
Subject  to 

teC 

jp  =  E  yt 

teC 

1  =  E  V 

teC 

po  “b  Mfhf  ^  dj^s  X  dt^r 

0<Xf 

hte  {0,1} 


Vjo  G  P, 

'ip  G  P, 

ip  G  P, 

Vt  G  T, 
it  G  G  P, 
Vt  G  P. 


(2.5a) 

(2.5b) 

(2.5c) 

(2.5d) 


The  model  forces  the  positions  of  the  source  and  receiver  to  be  inside  the  convex  hull 
by  defining  them  as  linear  combinations  of  the  vertices  in  Equations  (2.5a)  to  (2.5c).  We 
introduce  the  binary  decision  variable  ht  for  each  target  t  G  P  specifying  in  Equation  (2.5d) 
whether  a  target  t  is  detected,  ht  =  0,  or  hidden,  ht  =  1.  If  t  is  not  covered,  then  a  large 
number  Mt  is  added  to  the  range  of  the  day,  po,  making  sure  the  constraint  is  still  feasible. 

Camm,  Raturi,  and  Tsubakitani  (1990)  show  that  tight  penalty  parameters  like  Mt  reduce 
the  time  a  solver  needs  to  find  an  optimal  solution.  In  order  to  make  sure  the  INEP  is  always 
feasible,  the  left-hand  side  (EHS)  of  this  equation  has  to  be  equal  to  the  largest  product  of 
distances  possible  for  each  target.  Since  the  largest  distance  from  a  target  position  to  a 
point  in  a  convex  hull  is  to  one  of  its  vertices,  we  can  set  Mt  to  its  square.  We  can  also 
substract  po  since  it  is  already  added  on  the  EHS.  Thus  we  have 


Mt  =  m&xd}fi 

t’ec  ’ 


-Po- 


(2.6) 
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It  is  important  to  notice  that  there  is  generally  more  than  one  optimal  solution  with  this 
model.  First  of  all,  interehanging  the  souree  and  reeeiver  in  a  bistatie  environment  does 
not  change  the  outcome.  The  Cassini  oval  drawn  by  both  sensors  is  still  the  same.  But 
more  importantly,  if  a  sensor  is  moved  slightly  in  one  direetion,  it  may  still  deteet  the  same 
targets.  And  last  but  not  least,  there  might  be  a  solution  deteeting  different  targets,  but  the 
same  number  of  targets.  Henee,  in  most  seenarios  there  exist  infinitely  many  solutions. 

2.3.2  Other  Networks 

The  model  deseribed  in  Seetion  2.3.1  only  supports  bistatie  eookie  eutter  sensor  networks. 
Consequently,  there  are  two  ways  to  extend  this  model:  finding  positions  for  multiple 
sourees  and  reeeivers,  and  expanding  to  diffuse  sensor  models.  Though  models  ean  be 
formulated,  up  to  this  point  we  are  not  able  to  define  a  model  that  runs  in  reasonable  time. 
In  a  future  study,  advaneed  versions  of  this  model  ean  be  further  investigated. 
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CHAPTER  3: 

Divide  Best  Sector  Algorithm 


In  this  chapter  we  slightly  change  the  scenario.  In  an  operational  setting,  it  is  not  unlikely 
that  receivers  are  already  deployed  throughout  a  field  by  plane  or  ship.  As  soon  as  a  hostile 
underwater  asset  is  suspected  in  the  area  of  responsibility  (AoR),  e.g.,  by  passive  detection, 
sources  are  used  to  send  out  pings  in  order  to  find  the  exact  position  of  the  intruder  (Navy 
personnel,  personal  communication,  April  2014).  As  discussed,  it  is  crucial  to  send  out  as 
few  sources  and  subsequent  pings  as  possible. 

3.1  Problem  Statement 

Consider  the  following  scenario:  there  is  a  set  T  of  targets  with  fixed  positions  that  have 
to  be  monitored  or  detected.  Furthermore,  there  is  a  set  R  of  receivers  whose  positions  we 
already  fixed.  The  objective  is  to  find  optimal  positions  for  a  limited  number  of  sources, 
defined  in  set  S.  We  do  not  assume  a  particular  sensor  model;  rather,  we  discuss  a  method 
that  works  with  all  sensor  models  shown  in  Figure  1.3. 

Figure  3.1  shows  a  generic  scenario  with  fixed  targets  and  receivers.  The  z-axis  represents 
the  average  detection  probability  assuming  a  single  source  is  deployed  at  the  respective 
position.  Multiple  local  maxima  in  this  plot  show  that  the  objective  function  is  nonconvex. 
Nonetheless  we  now  develop  an  algorithm  guaranteed  to  solve  this  problem  to  near  global 
optimality  when  |5|  =  1. 

3.2  Algorithm  Development 

In  this  section  we  develop  the  Divide  Best  Sector  (DiBS)  algorithm.  The  algorithm  par¬ 
titions  the  area  of  possible  solutions  into  sectors.  At  each  iteration  the  sector  with  the 
highest  upper  bound  for  the  objective  function  is  further  divided  into  smaller  sectors.  In 
Section  2.1  we  showed  that  optimal  sensor  positions  can  always  be  found  inside  the  convex 
hull  spanned  by  the  targets.  We  slightly  relax  this  condition  to  the  smallest  rectangle  with 
edges  parallel  to  the  v:-  and  y-axes  that  contains  all  targets.  This  way  we  can  easily  divide 
the  area  of  possible  source  locations  into  smaller  sectors  by  slicing  horizontally  as  well  as 
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Figure  3.1:  Problem  Statement  -  The  optimal  source  position  with  fixed  target  and  receiver 
locations  needs  to  be  found.  The  z-axis  represents  the  average  detection  probability  assuming  a 
source  is  deployed  at  the  respective  position.  Multiple  local  maxima  indicate  nonconvexity. 
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vertically. 

Each  sector  y  of  the  resulting  set  of  sectors  E  now  is  evaluated  with  respect  to  an  up¬ 
per  bound,  uB(7),  for  the  objective  function.  This  is  the  heart  of  DiBS  and  illustrated  in 
Figure  3.2.  Since  receivers  and  targets  occupy  fixed  positions  in  this  scenario,  the  target- 
receiver  distance  di^r  is  constant.  Thus,  the  equivalent  range  pt,s,r  and  subsequently  the 
detection  probability  Pt  only  depend  on  the  target-source  distance  dt^s-  Hence,  for  every 
target  we  determine  the  hypothetical  source  position  that  is  closest  to  the  target  but  still 
inside  y.  These  hypothetical  source  positions  are  either  on  the  edge  of  the  sector  or  the 
target  position  itself,  as  seen  in  Figure  3.2. 

With  the  shortest  distance  dt^j  from  target  t  to  its  respective  hypothetical  source  position 
inside  sector  y  we  can  compute  the  highest  possible  detection  probability  Pt,  assuming  the 
source  location  is  inside  y.  Having  done  this  with  all  targets,  we  compute  the  upper  bound, 
uB  (7) ,  according  to  Table  1.1.  The  following  pseudocode  finds  the  upper  bound  for  a  given 
sector  7,  where  x™"  (y™ '")  is  the  lowest  and  the  highest  x  (y)  coordinate  for  7. 
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Figure  3.2:  Determine  Upper  Bound  -  For  each  target  t  there  exists  a  hypothetical  source 
position  inside  sector  /  that  has  the  shortest  distance  dtj  to  the  target.  This  distance  determines 
the  highest  possible  detection  probability  with  a  source  position  inside  y.  Merging  the  probabilities 
for  all  targets  results  in  an  upper  bound  for  y. 


1:  procedure  Upper. Bound(7) 

2:  for  allt  eT  do 

3 :  ^  min  (max  (x™" ,  Xf ) , 

4:  ^  min(max(3;™'',yf),y™“) 

5:  dt^Y  ^  \/{^t-XsY  +  {yt-ysY 

6:  compute  Pt 

7:  end  for 

8:  compute  objective  value  Z 

9:  return  Z 

10:  end  procedure 


>  use  chosen  sensor  model 
>  use  chosen  objective 


In  a  next  step  we  pick  the  sector  7  G  F  in  which  uB(7)  is  maximized  and  divide  it  into 
smaller  sectors  with  respective  upper  bounds.  Consistently  repeating  this  method  results  in 
smaller  sectors  and  subsequently  tighter  upper  bounds  until  a  termination  condition  is  met. 
Various  termination  conditions  are  conceivable,  such  as  a  maximum  sector  size,  optimality 
range,  etc.  The  following  pseudocode  shows  the  workings  of  DiBS. 
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1;  procedure  DiBS 
2:  create  initial  F 

3:  optimal  ^  FALSE 

4:  while  -^optimal  do 

5;  select  7  G  r  such  that  Upper.Bound(7)  is  maximal 

6:  if  termination  condition  is  met  then  >  use  chosen  termination  condition 

7;  optimal  ^  TRUE 

8:  else 

9:  r'  ^  split  7  into  smaller  sectors 

10:  r^rur\{7} 

11:  end  if 

12:  end  while 

13:  return  7 

14:  end  procedure 

Eigure  3.3  shows  the  example  from  Eigure  3.1  solved  with  DiBS.  It  vividly  illustrates  how 
sectors  with  a  low  upper  bound  are  ignored  while  promising  areas  are  further  explored.  Eor 
the  rest  of  the  chapter  we  look  at  various  details  of  DiBS  and  give  recommendations  on 
how  to  use  it  most  efficiently. 


3.3  Termination  Conditions 

Eirst,  we  want  to  make  sure,  is  that  the  algorithm  eventually  finishes  when  it  meets  a  specific 
condition.  This  section  discusses  various  termination  conditions  based  on  sector  size  and 
optimality  gap. 

3.3.1  Termination  By  Sector  Size 

There  are  two  ways  to  define  the  sector  size:  by  its  area  and  by  its  edge  lengths.  A  small 
area,  however,  can  be  the  result  of  a  very  narrow  but  long  rectangle.  This  metric  is  not  very 
useful,  because  detection  probabilities  can  change  greatly  along  a  line,  as  well  as  in  a  long 
rectangle.  Thus,  even  if  a  sector  with  small  area  has  a  high  upper  bound,  a  sensor  placed  in 
this  sector  may  perform  badly.  Additionally,  telling  a  decision  maker  to  place  the  source  in 
a  0.01  X  10, 000  foot  rectangle  might  not  be  useful  in  practice. 
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Figure  3.3:  Algorithm  Example  -  A  scenario  with  ten  targets  and  three  receivers  is  given  (see 
Figure  3.1).  DiBS  finds  the  source  position  with  the  highest  average  detection  probability  by 
dividing  only  promising  sectors,  i.e.,  sectors  with  a  high  upper  bound,  into  smaller  ones. 


On  the  other  hand,  using  a  maximum  edge  length  strongly  eonstrains  possible  souree  loea- 
tions.  Choosing  a  maximum  edge  length  of  ten  feet,  for  example,  results  in  a  final  seetor 
that  has  size  of  at  most  10x10  feet,  a  suffieient  preeision  for  helieopter  dips.  Moreover  the 
ehange  in  deteetion  probabilities  in  sueh  a  region  is  very  limited. 

Therefore,  we  disregard  using  area  as  the  termination  eondition.  For  the  rest  of  this  ehapter, 
termination  by  seetor  size  denotes  termination  when  the  length  of  the  longest  edge  of  the 
seetor  with  the  highest  upper  bound  is  suffieiently  small. 

3.3.2  Termination  By  Optimality  Gap 

Though  simple  and  useful,  it  turns  out  that  termination  by  seetor  size  has  a  drawbaek.  If 
we  request  a  seetor  with  very  short  edges,  the  differenee  for  the  upper  bounds  between 
adjaeent,  same-sized  seetors  beeomes  negligible.  This  leads  to  dividing  all  seetors  with  the 
same  size  in  an  area  before  eontinuing  with  smaller  ones.  This,  however,  takes  the  majority 
of  the  time  and,  on  top  of  this,  the  gain  in  objeetive  value  is  almost  always  negligible. 

Henee,  we  expand  the  termination  eondition  with  a  teehnique  used  by  solvers:  we  introduee 
a  tolerated  optimality  gap,  £.  By  allowing  DiBS  to  stop  as  soon  as  it  finds  a  speeifie  sensor 
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position  s'  that  has  an  objective  value  within  e  of  the  upper  bound,  we  are  able  to  stem 
the  described  effect.  The  effectiveness  of  the  optimality  gap  strongly  depends  on  which 
position  we  choose  for  s'. 

We  take  a  closer  look  at  the  methods  summarized  by  Table  3.1.  A  method  with  fewer 
potential  source  positions,  like  Center,  does  not  require  high  computational  power  but  is 
subject  to  missing  the  optimality  gap  more  often.  On  the  other  hand,  checking  five  positions 
with  Corners+ Center  is  more  precise  but  increases  run  time.  The  last  method.  Targets,  is 
an  add-on  for  the  previous  ones.  The  idea  is  that  the  detection  probability  for  a  single  target 
is  replacedmaximizedmaximum  only  if  at  least  one  sensor  is  at  the  target’s  position.  Since 
there  are  sectors  that  do  not  contain  any  targets,  it  can  only  be  used  in  conjunction  with 
another  method. 


Method 

Amount 

Remarks 

Center 

1 

uses  seetor’s  center  point 

Corners 

4 

uses  sector’s  vertices 

Corners-i-Center 

5 

combines  Corners  and  Center 

Targets 

0,1,... 

uses  target  positions  inside  sector 

Table  3.1:  Optimality  Gap  Methods  -  These  Methods  look  for  a  specific  source  position  inside 
a  sector  in  order  to  determine  the  optimality  gap.  The  last  method,  Targets,  has  to  be  used  in 
conjunction  with  one  of  the  other  methods  since  sectors  with  no  targets  inside  exists. 


The  example  in  Figure  3.4  shows  how  an  optimality  gap  limits  the  number  of  created  and 
evaluated  sectors  and  thus  the  run  time.  In  this  particular  problem  instance  the  baseline 
case  created  more  than  16,000  sectors  for  a  demanded  precision  of  0.001.  Utilizing  an 
optimality  gap  of  £  =  0.05,  however,  already  terminates  DiBS  after  about  200  created 
sectors  with  only  a  minor  difference  between  Center  and  Comers. 

Running  more  examples,  we  observe  that  the  gap  between  Center  and  Comers  always  is 
small.  Most  of  the  time  Comers  terminates  with  fewer  sectors  than  Center,  because  it  can 
choose  from  four  potential  source  positions  instead  of  just  one.  Center,  by  contrast,  has 
the  better  run  time.  Since  Corners-i-Center  combines  both  methods  it  always  uses  the  same 
number  of  sectors  as  the  method  with  fewer  sectors,  but  it  takes  the  longest  run  time.  In 
all  tests  we  ran,  we  never  observed  that  adding  the  method  Targets  had  any  effect  on  the 
outcome  except  for  increasing  run  time. 
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Figure  3.4:  Optimality  Gap  Example  -  While  the  number  of  evaluated  sectors  for  the  baseline 
case  without  optimality  gap  feature  increases  with  higher  demanded  precision,  using  an  optimality 
gap  terminates  the  algorithm  earlier. 

Including  the  termination  eondition  with  an  optimality  gap  is  strongly  reeommended  sinee 
it  signifieantly  reduees  the  number  of  created  and  evaluated  sectors,  especially  for  high 
precision  calculations.  Even  though  the  run  time  per  iteration  slightly  inereases,  the  total 
run  time  ean  reduee  to  a  fraetion  of  the  original  run  time,  whereas  the  ehoiee  of  method  for 
finding  the  definite  souree  position  s'  hardly  ehanges  the  outeome. 

3.3.3  Utilizing  Lower  Bounds 

On  eloser  examination,  we  realize  that  the  methods  described  in  the  previous  section  pro¬ 
vide  a  lower  bound,  IB (7).  This  section  discusses  whether  we  ean  use  lower  bounds  to 
disregard  seetors.  Additionally,  we  show  the  performanee  of  DiBS  using  lower  bounds. 

Consider  disearding  a  hypothetieal  seetor  /  whose  upper  bound,  uB(/),  is  less  or  equal 
than  the  lower  bound,  IB (7),  of  another  seetor  7.  Surprisingly,  this  is  not  neeessary.  Sinee 
IB  (7)  is  the  objective  value  for  a  definite  souree  position  in  7,  it  is  always  included  in  one 
seetor  of  7’s  split.  Henee  there  is  always  a  seetor  out  of  7’s  split  that  still  has  a  higher  upper 
and  lower  bound  than  /.  Thus,  /  will  never  be  chosen  as  the  best  seetor,  and  therefore, 
disearding  does  not  reduee  the  total  number  of  sectors  created. 

It  does,  however,  reduce  the  number  of  sectors  in  a  list  that  has  to  be  ordered  by  upper 
bounds.  Most  programming  languages  use  fast  sorting  algorithms  like  Quicksort.  So  or- 


29 


dering  a  large  list  eompared  to  other  operations  in  DiBS  is  negligible. 

Nevertheless,  we  ean  use  the  lower  bound  to  indireetly  measure  the  goodness  of  the  upper 
bound.  A  tight  upper  bound  is  important  for  faster  identifieation  of  a  solution.  Figure  3.5 
illustrates  the  change  of  gap  between  upper  and  lower  bound  throughout  the  execution  of 
DiBS.  Here,  the  right  plot  shows  the  underlying  trend  of  a  decreasing  gap  the  longer  the 
algorithm  runs.  Jumps  to  lower  values  for  the  lower  bound  occur  when  the  next  sector 
with  highest  upper  bound  is  much  larger  than  the  previous  sector  and  therefore  can  provide 
a  bigger  gap  between  lower  and  upper  bound.  The  left  plot  verifies  our  suspicion  that 
smaller  sectors  yield  smaller  gaps  and  therefore  tighter  upper  bounds. 
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Figure  3.5:  DiBS  Goodness  -  With  increasing  number  of  iterations  the  gap  between  upper  and 
lower  bound  decreases  (left  plot).  Steep  cuts  on  lower  bound  values  occur  on  iterations  where  a 
small  sector  is  sufficiently  split  and  DiBS  jumps  back  to  a  larger  sector  and  thus  larger  gap.  The 
right  plot  mirrors  the  relationship  between  longest  edge  size  as  a  multiple  of  longest  initial  edge 
and  gap  size.  Smaller  sectors  yield  tighter  upper  bounds. 


Summarizing,  we  conclude  that  lower  bounds  do  not  accelerate  DiBS.  They  can,  however, 
be  used  to  show  that  upper  bounds  are  getting  tighter  along  the  algorithm.  Hence,  the 
algorithm  quickly  achieves  a  small  optimality  gap,  which  can  be  used  to  terminate  DiBS 
early  as  shown  in  Section  3.3.2. 


3.4  Sector  Splitting 

Since  the  algorithm  begins  with  creating  the  initial  set  of  sectors  F,  this  procedure  and 
subsequently  the  dividing  into  smaller  sectors  is  the  next  feature  of  DiBS  into  which  we 
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take  a  closer  look.  The  question  is  how  many  vertical  and  horizontal  slices  are  efficient 
with  respect  to  complexity  and  run  time?  First  we  investigate  a  fixed  number  of  slices  and 
expand  later  on  to  dynamic  splitting  methods. 

3.4.1  Fixed  Splitting 

More  slices  reduce  sector  sizes  faster  and  therefore  meet  the  termination  condition  in  fewer 
iterations  of  DiBS.  Each  iteration,  however,  requires  more  computations  and  thus  more 
time.  The  plots  in  Figure  3.6  show  the  algorithm’s  performance  for  a  particular  problem 
instance  using  various  settings  for  initial  division  and  further  splitting.  The  derived  source 
position  is  always  the  same  for  all  settings.  As  assumed,  the  number  of  iterations  decreases 
quickly  with  more  initial  sectors  and  even  more  with  a  higher  number  of  new  sectors  per 
iteration.  The  positive  effect,  however,  stagnates  very  quickly. 

number  of  initial  sectors: - 4; - 36;  100;  225; - 400 


new  sectors  per  iteration  new  sectors  per  iteration 

Figure  3.6:  Sector  Splitting  -  Algorithm’s  performance  with  different  splitting  settings.  Splitting 
into  many  sectors  per  iteration  increases  the  total  number  of  sectors  and  subsequently  the  run 
time  of  the  algorithm.  The  size  of  the  initial  set  of  sectors  also  affects  run  time  especially  for 
high  numbers  of  new  sectors  per  iteration.  The  reduction  in  number  of  total  iterations  is  not  as 
decisive. 

On  the  other  hand,  the  right  plot  in  Figure  3.6  reflects  the  increasing  number  of  created 
sectors.  While  the  change  in  the  number  of  iterations  is  between  8  and  62,  the  values  for 
the  number  of  created  sectors  range  from  270  to  over  10,000.  Running  DiBS  on  a  test 
machine  shows  that  the  run  time  is  directly  related  to  the  number  of  created  sectors,  e.g., 
0.9  seconds  for  an  initial  and  further  split  of  4  sectors  versus  28.7  seconds  for  an  initial  and 
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further  split  of  400  sectors.  The  number  of  iterations  does  not  seem  to  have  any  effect  on 
the  run  time. 

The  results  are  decisive  for  newly  created  sectors:  dividing  the  best  sector  into  fewer  new 
sectors  greatly  decreases  run  time  and  improves  algorithm  efficiency.  A  reason  for  this 
observation  lies  in  dividing  small  sectors.  We  know  from  Section  3.3.3  that  small  sectors 
have  tighter  upper  bounds.  Hence,  the  variance  in  upper  bounds  of  their  split  sectors  is 
small.  So  after  splitting  a  small  sector  into  a  set  of  many  small  sectors  F',  it  is  very  likely 
that  the  algorithm  has  to  evaluate  and,  even  worse,  split  most  of  the  sectors  in  F'  before  it 
is  able  to  continue  with  a  smaller  sized  sector.  This  behavior  creates  many  sectors  that  are 
unnecessarily  evaluated  and  therefore  waste  resources.  A  different  approach  is  discussed 
in  Section  3.4.2. 

The  creation  of  the  initial  set  of  sectors  also  has  an  effect  on  the  total  number  of  created 
sectors.  Especially  at  high  numbers  of  new  sectors  per  iteration  it  becomes  a  critical  factor, 
as  seen  in  Figure  3.6.  It  appears  a  smaller  initial  set  performs  more  efficiently. 

In  most  cases,  starting  with  a  small  set  of  sectors  and  slicing  the  best  sector  into  a  small 
number  of  new  sectors  is  most  efficient  and  performs  best.  We  recommend  to  always  start 
with  this  setting  and  expand  as  circumstances  require. 

3.4.2  Dynamic  Splitting 

We  discovered  in  Section  3.4.1  that  splitting  into  many  sectors  per  iteration  produces  a  lot 
of  overhead,  especially  at  small  sector  sizes.  The  question  arises  whether  one  can  reduce 
the  number  of  created  sectors  via  dynamic  splitting.  In  dynamic  splitting,  large  sectors  or 
sectors  with  a  high  optimality  gap  are  split  into  more  sectors  than  small  sectors  or  sectors 
with  tighter  upper  bounds.  We  leave  that  question  open  for  future  work. 

3.4.3  Plane  Rotation 

The  DiBS  algorithm  starts  with  finding  the  smallest  rectangle  with  edges  parallel  to  the 
X-  and  y-axes  that  contains  all  targets.  We  know  from  Section  2.1  that  the  optimal  sensor 
position  lies  inside  the  convex  hull  spanned  by  the  targets.  Hence,  it  is  possible  that  the 
resulting  rectangle  contains  a  large  overhang  which  is  subject  to  being  evaluated  unneces¬ 
sarily.  This  section  discusses  the  use  of  the  minimum-area  rectangle  without  constrained 
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edges  and  how  to  rotate  the  plane,  such  that  DiBS  still  can  be  applied. 

Freeman  and  Shapira  (1975)  prove  that  the  minumum-area  rectangle  encasing  a  convex  hull 
has  a  side  collinear  with  one  of  the  edges  of  the  convex  hull,  like  the  examples  in  Figure  3.7. 
Hence  we  choose  the  rectangle  with  smallest  area  out  of  the  |C|,  i.e.,  the  number  of  vertices 
of  the  convex  hull,  rectangles  we  can  form.  Using  the  Graham  scan  to  find  C  has  the 
advantage  that  the  vertices  are  already  ordered  counter  clockwise.  Therefore,  we  can  use 
this  output  to  easily  construct  the  edges  of  the  convex  hull. 


Figure  3.7:  Minimum-Area  Rectangle  Example  -  The  plot  shows  the  original  rectangle  as 
well  as  two  rectangles  that  have  one  side  collinear  with  one  side  of  the  convex  hull.  The  Tg  — Tg 
rectangle  is  the  minimum-area  rectangle  that  encases  all  targets. 

The  example  in  Figure  3.7  shows  the  original  rectangle  produced  by  DiBS  as  well  as  the 
rectangles  with  the  smallest  and  the  largest  area  that  have  one  side  collinear  with  one  of  the 
edges  of  the  convex  hull.  The  respective  areas  in  square  units  are 


original:  636,  T(,  -  T^:  505.375,  Tio  -  Ti:  619.420. 


This  example  shows  that  using  the  minimum-area  rectangle  can  greatly  reduce  the  over¬ 
hang.  Having  determined  the  minimum-area  rectangle,  we  need  to  rotate  the  plane  such 
that  the  edges  of  the  rectangle  are  parallel  to  x-  and  y-axis.  An  efficient  way  to  do  this  is 
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applying  a  linear  transformation  to  the  sets  of  target  and  reeeiver  locations  (Leon,  2010, 
p.  166ff.).  Let  to  —  h  be  the  edge  that  is  collinear  with  the  minimum-area  rectangle.  Then 
the  linear  transformation  is  represented  by  the  matrix 


A  = 


cosO 
—  sin0 


sin0 

COS0 


where  9  =  tan  ^ 


yti 

Xt,  -XtJ 


The  matrix  multiplication  AX  provides  the  new  coordinates  for  targets  and  receivers,  where 
the  first  and  second  row  of  X  holds  the  x  and  y  coordinates  respectively. 


Since  plane  rotation  can  have  a  huge  impact  on  the  size  of  the  overhang,  it  should  always 
be  considered  before  running  DiBS.  There  is  also  the  fact  that  Graham  scan  and  linear 
transformation  are  very  efficient  and  therefore  executed  quickly. 


3.4.4  Handling  Uneven  Rectangles 

The  next  improvement  on  sector  splitting  aims  for  rectangles  with  uneven  edge  lengths. 
Obviously  the  short  edge  reaches  its  termination  condition  faster  than  the  long  edge  if  DiBS 
performs  the  same  number  of  slices  in  both  directions  at  each  iteration.  Since  termination 
by  sector  size  considers  only  the  longest  edge,  the  algorithm  continues  dividing  sectors  in 
both  directions.  Hence  it  creates  more  sectors  than  needed. 

There  are  two  ways  to  address  this  problem:  stopping  splits  in  directions  that  meet  the  ter¬ 
mination  condition  and  finding  an  initial  split  that  has  only  square  sectors.  The  first  method 
is  easy  to  implement  by  setting  the  number  of  slices  for  edges  that  meet  the  termination 
condition  to  zero.  There  is,  however,  already  a  problem  with  uneven  rectangles  before  the 
short  edge  meets  the  termination  condition.  As  opposed  to  squares  with  the  same  area, 
especially  long  rectangles  are  more  likely  to  have  a  big  variance  in  detection  probabilities 
(see  Section  3.3.1).  This  subsequently  leads  to  looser  upper  bounds. 

The  second  approach  bypasses  this  issue  by  starting  with  squares  at  the  first  iteration.  We 
already  discovered  in  Section  3.4.1,  that  a  small  set  of  initial  sectors  performs  more  effi¬ 
ciently.  So  we  are  looking  for  a  method  that  creates  a  low  number  of  squares  that  covers  an 
uneven  rectangle. 

Figure  3.8  shows  two  methods  to  create  initial  squares.  For  the  rest  of  this  section,  we 
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Figure  3.8:  Creating  initial  squares  -  The  first  method  uses  half  of  the  long  edge  b  for  square 
edge  length,  while  method  2  uses  the  short  edge.  The  red  shaded  area  reflects  the  created 
overhang  when  using  only  squares  for  the  initial  split. 


assume  that  a  and  b  are  the  short  and  the  long  edge  respeetively.  The  first  method  outs 
the  long  side  b  in  half  and  takes  the  result  as  the  edge  length  for  the  squares.  The  seoond 
method,  on  the  other  hand,  uses  the  short  edge  a  as  the  edge  length  for  the  squares.  Both 
methods,  however,  oreate  an  overhang,  shaded  red  in  the  figure,  that  hold  non-optimal 
solutions.  The  overhangs  are  oaloulated  as  follows: 

For  method  1:  —  aZ?.  (3.1) 

fO  if  a  I  h, 

For  method  2:  0.2=  \  (3.2) 

[[a  —  {b  mod  a)]a  otherwise. 

In  order  to  keep  the  unneeded  area  minimal,  we  always  ohoose  the  method  with  the  smallest 
overhang.  Table  3.2  lists  two  examples  with  different  ehoiees  of  the  method  used. 


a 

b 

r  =  ajb 

Method  1  Overhang 

Method  2  Overhang 

Select 

4 

5 

0.8 

52-4x5  =  5 

[4  -  (5  mod  4)]4  =  12 

Method  1 

2 

5 

0.4 

52-2x5  =  15 

[2  -  (5  mod  2)]2  =  2 

Method  2 

Table  3.2:  Examples  for  Overhang  -  For  the  first  example,  method  1  has  the  smaller  overhang, 
hence  we  select  it.  For  the  second  example,  however,  we  choose  method  2  because  of  its  smaller 
overhang. 


Going  a  step  farther,  we  want  to  see  whether  we  are  able  to  decide  which  method  to  use  by 
solely  looking  at  the  ratio  between  long  and  short  edge. 
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Theorem  3.1.  Method  1  has  less  overhang  than  method  2  if  and  only  if  ^<r  <\,  where 
r  =  alb. 


Proof  Without  loss  of  generality,  we  assume  b  =1.  Then  0  <  r  =  a  <  1 .  Furthermore  the 
overhangs  for  both  methods  compute  as  follows: 

0.1  =  b^  —  ab  =  I  —  r^ 

fO  if  alb, 

^2={ 

[[a  —  {b  mod  a)]a  =  [r  —  (1  mod  r)]r  otherwise. 

Method  1  has  less  overhang  if  and  only  if  <  0,  which  is  visualized  in  Figure  3.9. 

With  a  ratio  of  r  =  1,  i.e.,  both  sides  have  the  same  length,  both  methods  do  not  have  an 
overhang  since  a  and  b  already  form  a  square.  Discontinuities  occur  on  positions  where  a 
divides  b,  such  as  Yi,  V-i,  14  and  so  on.  Here,  method  2  produces  no  overhang.  Moreover,  by 
inspection,  we  see  —  ^^2  >  0  for  0  <  r  <  0.5. 


Figure  3.9:  Theorem  3.1  Proof  -  The  figure  shows  the  difference  between  method  I's  overhang, 
Hi,  and  method  2's  overhang,  0.2,  with  respect  to  the  ratio  r  =  alb.  A  positive  value  indicates 
a  smaller  overhang  for  method  2,  while  method  1  is  preferred  at  negative  values. 
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Hence,  we  are  looking  for  a  root  in  0.5  <  r  <  1,  where  a\b  and  1  mod  r  =  1  —  r.  Thus 


—  ^^2  =  1  —  r  —  [r  —  (1  mod  r)]r 
=  l-r-[r-(l-r)]r 
=  l-2r2, 

of  which  the  only  nonnegative  root  is  TTT  □ 


When  DiBS  splits  a  sector,  it  does  not  recognize  whether  or  not  a  newly  created  sector  is 
in  the  overhang.  An  additional  check  removes  those  sectors  and  prevents  selecting  them  as 
best  sector  for  further  splitting. 


3.5  Finding  Multiple  Sensor  Locations 

A  way  to  increase  the  number  of  possible  applications  of  DiBS  is  increasing  the  number 
of  sensor  positions  optimized.  This  section  investigates  an  iterative  method  as  well  as  a 
simultaneous  approach. 


3.5.1  Placing  Sources  Iteratively 

We  consider  a  set  S  of  sensors,  that  has  to  be  optimally  deployed  in  a  scenario  with  fixed 
targets  and  receivers.  Furthermore,  we  assume  a  cookie  cutter  sensor  model.  The  position 
for  the  first  source  is  determined  by  run  DiBS  without  modifications.  In  a  next  step  we 
remove  all  targets  that  are  detected  by  the  algorithm’s  solution  and  running  DiBS  again  for 
the  second  sensor  position.  We  repeat  those  step  until  either  all  targets  are  detected  or  all 
sources  are  set.  The  following  pseudocode  illustrates  the  iterative  method. 
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1;  procedure  Iterative. Method 
2:  while  |r|  >  0  A  |5|  >  0  do 

3:  select  s  e  S 

4:  run  DiBS 

5;  T  ■^T\  {detected  targets} 

6:  5^5\{4 

7:  end  while 

8;  end  procedure 

Since  we  place  sources  in  a  greedy  approach,  the  resulting  solution  generally  is  not  optimal. 
Notwithstanding,  this  algorithm  reveals  some  information  about  the  optimal  solution.  First 
of  all,  in  case  all  targets  are  detected  by  this  method,  the  solution  indeed  is  optimal. 

Otherwise,  we  can  use  the  outcome  to  define  a  lower  and  upper  bound  for  this  problem 
by  applying  an  observation  made  by  Craparo,  How,  and  Modiano  (2011).  Leveraging  sub¬ 
modularity  in  a  similar  coverage  setting,  the  authors  determine  that  an  iterative  algorithm 
has  an  optimality  gap  of  at  most  1  —  e~^  ~  0.632.  Hence,  we  define  the  lower  bound,  IB, 
as  the  outcome  of  the  iterative  approach  and  the  upper  bound 


^  [1.582  xlBj. 

We  conclude  that  even  though  the  iterative  method  does  not  guarantee  an  optimal  solution, 
it  is  a  quick  way  to  get  an  estimate  about  the  quality  of  other  solutions.  Combined  with  our 
observations  from  Section  2.2,  we  are  able  to  narrow  down  the  optimal  solution  even  more. 

3.5.2  Placing  Sources  Simultaneously 

In  contrast  to  the  previous  method,  a  different  approach  is  placing  multiple  sources  simul¬ 
taneously.  In  doing  so,  we  place  each  source  in  an  arbitrary  sector  and  evaluate  upper 
bounds  of  occupied  sectors.  Merging  those  upper  bounds  yields  the  value  for  the  objective 
function.  In  order  to  exhaust  all  possible  solutions,  we  evaluate  all  combinations  of  sec¬ 
tors.  We  have  to  take  into  account  that  multiple  sources  might  occupy  the  same  sector  and 
that  the  order  of  sources  is  indifferent.  Hence,  by  definition  the  number  of  combinations  is 
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calculated  by  combinations  of  multisets  (Brualdi,  2010,  p.  52ff.)  as 


/|s|  +  |r|-i\ 

(bi  +  in-i)! 

V  |5|  J 

|X|!(|r|-l)! 

where  |5|  is  the  number  of  sensors  and  |r|  the  number  of  sectors  in  a  particular  iteration. 
The  base  ease  with  one  souree  represents  the  original  DiBS  where  the  number  of  eombina- 
tions  equals  the  number  of  seetors,  i.e.,  plaeing  the  souree  onee  in  eaeh  seetor.  With  two 
sources  and  four  sectors  we  already  have  the  following  ten  eombinations. 

l:{7i,7i}  2:{7i,r2}  3:{7i,73}  4:{7i,74}  5:{72,72} 

6:{72,73}  7:{72,74}  8:{73,73}  9:{73,74}  10:{74,74} 

The  number  of  combinations  grows  quickly  with  the  number  of  sectors  and  sources  whieh 
is  visualized  in  Figure  3.10.  While  plaeing  just  a  few  sourees  simultaneously  might  still 
be  manageable,  using  this  approach  with  many  sensors  quiekly  becomes  an  issue.  This  is 
aggravated  by  the  fact  that  selecting  and  dividing  multiple  seetors  per  iteration  inereases 
the  number  of  seetors  and  subsequently  the  number  of  eombinations  even  more  quiekly. 
A  future  study  might  find  ways  to  reduce  that  number  or  find  another  method  to  deal  with 
multiple  sourees. 


number  of  sources: - 1; - 2; - 3; - 5; - 10 


Figure  3.10:  Sector  Combinations  -  This  plot  displays  the  number  of  combinations  that  have 
to  be  evaluated  with  respect  to  the  number  of  sectors  and  sources.  The  base  case  with  one 
source  represents  the  original  DiBS. 
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CHAPTER  4: 
Conclusion 


We  started  by  developing  some  insights  that  help  bound  the  problem.  We  showed  that  plac¬ 
ing  sensors  inside  the  convex  hull  that  encases  all  targets  results  in  detection  probabilities 
as  least  as  good  as  outside  the  hull.  This  important  observation  limits  the  area  where  we 
have  to  search  for  optimal  sensor  positions  regardless  of  the  sensor  model  we  use. 

Assuming  a  cookie  cutter  sensor  model,  we  can  utilize  clusters  of  RDCs.  The  number 
of  elements  in  the  minimal  set  of  clusters,  G,  represents  a  lower  bound  on  the  number  of 
sensors  required  to  cover  all  targets.  In  addition,  with  a  fixed  number  of  sensors,  n,  we 
are  able  to  determine  an  upper  bound  on  the  number  of  targets  by  selecting  n  clusters  from 
G,  such  that  the  number  of  contained  targets  is  maximized.  Other  circles,  such  as  RCs  and 
SCs,  are  means  to  narrow  down  possible  sensor  positions  assuming  the  location  of  one  kind 
of  sensor  is  fixed. 

In  a  subsequent  step,  we  formulated  an  INLP  to  find  optimal  positions  for  a  bistatic  cookie 
cutter  sensor  network.  Up  to  this  point,  expanding  this  model  to  multiple  sensors  or  other 
sensor  models  does  not  provide  a  solution  in  reasonable  time.  Hence,  we  developed  the 
DiBS  algorithm. 

DiBS  assumes  that  target  and  receiver  positions  are  fixed  and  determines  the  optimal  posi¬ 
tion  for  a  single  source  regardless  of  the  applied  sensor  model.  By  dividing  sectors  with  the 
highest  upper  bound  for  the  objective  function,  it  eventually  reaches  a  termination  condition 
like  maximum  size  of  a  sector’s  longest  edge  or  a  maximum  optimality  gap. 

Further  investigations  of  the  DiBS  algorithm’s  details  resulted  in  the  following  observations 
and  recommendations.  The  number  of  new  sectors  created  at  each  iteration  as  well  as  the 
number  of  initial  sectors  should  be  kept  small  in  order  to  reduce  the  total  number  of  created 
and  evaluated  sectors  and  subsequently  the  run  time  of  DiBS. 

Rotating  the  plane  such  that  the  edges  of  the  minimum-area  rectangle  encasing  all  targets 
are  parallel  to  the  x-  and  y-axis  results  in  a  smaller  needless  overhang.  Because  of  the  high 
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efficiency  to  find  the  minimum-area  rectangle  and  to  rotate  the  plane,  this  improvement 
should  always  be  considered. 

Furthermore,  we  discussed  two  ways  to  handle  uneven  rectangles  that  can  lead  to  creating 
and  evaluating  many  unnecessary  sectors.  Both  methods  approach  this  issue  differently. 
Their  efficiency  depends  on  the  given  scenario.  We  recommend,  however,  to  apply  one  of 
the  methods  to  reduce  the  algorithm’s  run  time. 

Finally,  we  illustrate  two  methods  to  apply  DiBS  to  problems  where  we  have  to  find  mul¬ 
tiple  source  locations.  The  faster  iterative  method  assuming  a  cookie  cutter  sensor  model 
provides  lower  and  upper  bounds  on  the  number  of  detected  targets.  The  simultaneous  ap¬ 
proach,  however,  quickly  creates  a  large  number  of  combinations  that  have  to  evaluated.  A 
future  study  can  develop  ways  to  reduce  the  number  of  combinations  or  find  other  means 
to  find  multiple  source  locations  using  DiBS. 
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APPENDIX  A: 

Excel  Analysis  And  Decision  Tool 


For  a  better  understanding  and  visualization  of  the  observations  from  Chapter  2  ,  we  devel¬ 
oped  an  analysis  and  decision  tool  in  Microsoft  Excel.  It  provides  a  number  of  features  for 
a  bistatic  cookie  cutter  sensor  model,  e.g.,  showing  the  Cassini  oval. 


The  Data  worksheet  shown  in  Figure  A.  1  contains  the  scenario  data  such  as  target  positions 
and  range  of  the  day.  Furthermore,  it  shows  the  positions  of  source  and  receiver  and  deter¬ 
mines  the  detection  probability  for  each  target  based  on  target-source  and  target-receiver 
distance. 


A 

B 

C 

D 

E  F  G  H  1  J  K  L 

1 

2 

Range  of  the  Day 

Object 

X 

V 

Src  Dist 

Rev  Dist 

detected 

3 

Po 

3.5 

T1 

-15 

-4.75 

18.6269294 

15.3121684 

FALSE 

4 

T2 

6.25 

-9.5 

13.9950884 

19.054199 

FALSE 

5 

Source 

X 

V 

T3 

-9.25 

-11 

18.2705364 

17.9669836 

FALSE 

6 

Sourcel 

1.6 

3.7 

T4 

0.75 

8.75 

5.12103505 

5.70306935 

FALSE 

7 

T5 

-4 

5.5 

5.88217647 

0.89442719 

TRUE 

8 

Receiver 

X 

V 

T6 

9 

-11.25 

16.681202 

22.0808175 

FALSE 

9 

Receiver  1 

-4.4 

6.3 

T7 

-9.5 

-15 

21.746264 

21.9020547 

FALSE 

10 

T8 

2.75 

11.5 

7.88431988 

8.84095583 

FALSE 

11 

UpdateTarget  Information 

1 

T9 

3.5 

-12.25 

16.0627675 

20.1621551 

FALSE 

12 

J 

TIO 

-4 

7.5 

6.76756973 

1.26491106 

TRUE 

13 

14 

LoadTargetstrom  .CSV 

1 

15 

16 

Figure  A.l:  MSN  Tool  Data  Sheet  -  The  Data  worksheet  contains  information  about  targets, 
sensors  and  range  of  the  day.  It  provides  methods  to  change  the  scenario. 


There  are  two  ways  to  change  the  target  data:  by  changing  the  values  in  the  table  or  by 
loading  a  comma-separated  values  (CSV)  file.  Using  the  first  method,  we  have  to  press 
the  Update  Target  Information  button  after  all  changes  were  made  in  order  to  apply  the 
changes.  The  second  method  initiated  by  clicking  Load  Targets  from  .CSV  does  not  require 
this  step.  The  following  example  shows  the  required  structure  for  a  CSV  file. 


11,-15,-4.75 

12,6.25,-9.5 


43 


Figure  A. 2  shows  the  Chart  worksheet  that  visualizes  the  scenario.  Here,  we  are  able  to 
show  or  hide  various  features.  The  scrollbars  on  the  top  and  on  the  right  move  either  source 
or  receiver  depending  on  what  is  chosen  on  the  left.  Detected  targets  are  represented  by 
empty  circles;  otherwise,  targets  are  represented  by  filled  circles.  The  Data  box  contains 
some  general  information  while  the  Selection  box  provides  a  short  report  about  each  ele¬ 
ment  in  the  chart  based  on  the  mouse  position. 


“  Options - 

[7  Show  Cassini  Oval 
W  Show  Convex  Hull 
W  Show  Range  of  the  Day  Circles 
r*  Show  Soufce  Circles 
^  Show  Receiver  Circles 
Use  Scrollbars  to 
O  niove  Source  ♦ 

(S  move  Receiver  A 

-  Data - 

Range  of  the  Day;  3.5 

Source  position:  (1.6;  3.7) 
Receiver  position:  (-4.4;  6.3) 

No.  of  covered  Targets;  2 

-  Selection - 

—  Range  of  the  Day  (T10) 

This  circle  intersects  with  the 
Range  of  the  Day  circle  of; 

T4.T5 

In  order  to  detect  T10^  at  least  one 
sensor  has  to  be  inside  this  circle. 


Find  Optimum 


4  ► 


Figure  A. 2:  MSN  Tool  Chart  -  The  scrollbars  move  either  source  or  receiver  based  on  the 
selection  on  the  right.  The  Selection  box  provides  additional  information  about  shown  objects  in 
the  chart  when  hovering  with  the  mouse.  In  this  example  target  Tio’s  RDC  is  selected. 

The  Find  Optimum  button  is  without  functionality  in  this  version  of  the  tool.  Here,  the 
INLP  from  Section  2.3  can  be  integrated  in  future  work. 
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APPENDIX  B: 

R  Implementation  Of  DIBS 


DiBS  <—  function  ( targets  ,  receivers  ,  rho_0  ,  longest_edge  ,  verbose  =  TRUE,  obj  =  ’  avg  ’  , 
sensor  =  ’  fermi  ’  ,  b  =  0.25,  gap_method  =  ’none’,  gap  =  0.05)  { 

#  Main  function  to  run  DiBS 

# 

i  krgs  : 

M  targets:  Data  frame  with  object,  x,  y  and  value. 

^  receivers :  Data  frame  with  object,  x  and  y. 

M  rho_0:  Range  of  the  day. 
i  long est _  edge :  Termination  condition. 

M  verbose:  If  TRUE,  prints  progress .  Default  is  TRUE. 

#  obj:  Objective.  Possible  values:  *avg’,  ^min\  Default  is  *avg\ 

i  sensor:  Sensor  model.  Possible  values:  *cookie*,  *fermi*,  ^ exp  ^ .  Default  is  *fermi*. 
M  b:  Diffusivity  parameter  for  Fermi  model.  Default  is  0.25. 

i  gap_method:  Optimality  gap  method.  Possible  values:  *none*,  ^center^,  *  corners  \  * 

corner  s  .  cent  er  * .  Default  is  ^none\ 

^  gap:  Accepted  optimality  gap.  Default  is  0.05. 

# 

#  Returns : 

#  Final  set  of  sectors,  where  the  first  row  yields  the  solution. 
terminate  <—  FALSE ;  iteration  <—  0 

data  <—  . create .  data ( targets  ,  receivers) 

sectors  <—  . create .  sectors (data  ,  rho_0  ,  obj,  sensor,  b,  ’initial’) 
while  (!  terminate  )  { 

iteration  <—  iteration  +  1 

sectors  <—  s e ctors [ order(~  s e c to r s $  upper . bound ,  sec  tors  $  longest . edge )  ,] 
gamma  <—  sectors  [1  ,] 

gamma$  lower  .  bound  <—  .  lower  .  bound  (  data  ,  gamma,  rho_0  ,  obj,  sensor,  b,  gap_method) 
gamma$  optimality  <—  gamma$ lower  .  bound  /  gamma$ upper  .  bound 
if  (gamma$ longe St  .  edge  <=  longest_edge  I  1  —  gamma$optimality  <=  gap) 
terminate  <-  TRUE 
else  { 

if  (verbose)  {cat(pasteO(’Iteration  ’,  iteration,  ':\n’));  print  (gamma)  ;  cat  (  ’  \  n  ’  )  } 
if  (gammaSxmax  —  gamma$xmin  <=  longest_edge )  xs  <—  c (gamma$xmin ,  gamma$xmax) 
else  xs  <—  seq  (gammaSxmin ,  gamma$xmax,  length,  out  =  3) 

if  (gammaSymax  —  gammaSymin  <=  longest_edge )  xs  <—  c (gamma$xmin ,  gamma$xmax) 
else  ys  <—  seq  (gammaSymin ,  gammaSymax,  length,  out  =  3) 

sectors  .prime  <—  .  create  .  sectors  (data  ,  rho_0  ,  obj,  sensor,  b,  xs  =  xs,  ys  =  ys, 

name  =  pasteO  (  ’  i  t  e  r  ’  ,  iteration)) 
sectors  <—  rbind  (  sectors  [  — 1  ,]  ,  sectors  .  prime  ) 


cat  ( pasteO  (’ Finished  after  ',  iteration  ,  ’  iterations  (evaluated  ’, 

nrow (sectors)  +  iteration  —  1,  ’  sectors)  :\n’)) 
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print  (gamma) 
return  ( sectors ) 


create  .  data  <—  function  ( targets  ,  receivers)  { 

^  Creates  the  data  frame  used  hy  DiBS  containing  target  locations , 

^  values  and  distances  to  receivers 

M  Args  : 

#  targets:  Data  frame  with  object,  x,  y  and  value. 

#  receivers :  Data  frame  with  object,  x  and  y. 

M  Returns : 

i  Created  data  frame. 

data  <—  t ar ge t s [, c (’ X y value ’) ] 

row .  names  (  data  )  <—  targetsSobject 

data[pasteO(’dist.’,  as.factor(receivers$object))]  <— 
sap  ply  (row  .  names  (  receivers),  function(r) 
supply  (row  .  names  ( targets),  function(t) 

sqrt(( targets [t ,  ’x’]— receivers [r ,  ’x’ ])^2+(targets [t ,  ’y’]— receivers [r ,  ’y’ ]) ^2))) 
return  ( data ) 


create  .  sectors  <—  function  ( data  ,  rho_0  ,  obj  ,  sensor  ,  b,  name,  xs  =  NUIX,  ys  =  NUIX)  { 

^  Creates  sectors  from  the  given  data  and  evaluates 
^  each  sector^s  longest  edge  and  upper  bound.  Uses  xs  and  ys  as 
i  slices.  If  xs  and  ys  are  NULL  then  it  creates  the  initial  split. 

# 

#  Args  : 

M  data:  Data  set  created  by  function  .  create .  data () . 
i  rho_0:  Range  of  the  day 

#  obj :  Objective. 

#  sensor:  Sensor  Model. 

M  b:  Diffusivity  parameter  for  Fermi  model, 

i  name:  identifier  for  the  set  of  sectors 
i  xs :  Vector  of  vertical  slices.  Default  is  NULL, 
i  ys :  Vector  of  horizontal  slices.  Default  is  NULL. 

# 

i  Returns : 

i  Data  frame  with  sector  coordinates  and  upper  bounds. 
if  ( is  .  null  (  xs  ) )  xs  <—  seq  ( min(  data$x)  ,  max(data$x),  length.  out=3) 
if  ( is  .  null  (  ys  ) )  ys  <—  seq  ( min(  data$y )  ,  max(data$y),  length.  out=3) 
sectors  <—  data.frame( 

sector  =  paste0(name,  1 :(( length  ( xs )  —  1)  *  (length(ys)  —  1))), 

xmin  =  rep  (  xs[— length  (  xs  )]  ,  length(ys)  —  1), 

xmax  =  rep(xs[  — 1],  length(ys)  —  1), 

ymin  =  rep  ( ys  [— length  (  ys  )]  ,  each  =  length(xs)  —  1), 

ymax  =  rep(ys[  — 1],  each  =  length(xs)  —  1)) 

s  e  c  tor  s  $  longe  s  t  .  edge  <—  pmax(  s  e  c  t  o  r  s  $xmax  —  sectors  $xmin  ,  sectorsSymax  —  s  ec  tor  s  $ymin ) 
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sectors  Supper  .  bound  <—  sapply  ( 1 :  nrow(  sectors  )  ,  function(s)  { 
switch(obj  ,  avg  =  value  <—  0,  min  =  value  <—  Inf, 

stop(pasteO(’objective  obj,  not  known  ’ ) ) ) 

if  (distr  ==  ’cookie’)  value  <—  0 
for(t  in  1 :  nrow(  data  ) )  { 

X  <—  min(max(  sectors  $xmin  [s],  data$x[t]),  sectors  $xmax  [  s  ] ) 
y  <—  min(max(  sectors  Symin  [s],  data$y[t]),  sectors  Symax  [  s  ] ) 

Pt  <—  .Pt(data,  t,  X,  y,  rho_0  ,  sensor,  b) 
if(sensor  ==  ’cookie’)  { 

value  <—  value  +  Pt*data[t,  ’value’] 

}  else  { 

switch  ( obj  , 

avg  =  value  <—  value  +  data[t,  ’value’]  *  Pt  /  nrow(data), 
min  =  value  <—  min(value  ,  data[t,  ’value’]  *  Pt) 

) 


return  (value  ) 

1) 

return  ( sectors ) 


.  Pt  <—  function  ( data  ,  t,  x,  y,  rho_0  ,  sensor,  b)  { 

i  Calculates  detection  probability  Pt  for  target  t  with  source  at  position  (x,  y) 

M  Args  : 

i  data:  Data  set  created  by  function  .  create .  data () . 

#  t :  target 

i  X :  source  *  s  x  coordinate 

#  y :  source  ^ s  y  coordinate 

#  rho_0:  Range  of  the  day 

M  sensor:  Sensor  Model. 

#  b:  Diffusivity  parameter  for  Fermi  model. 

# 

#  Returns : 

#  Detection  Probability  Pt 

d.ts  <—  sqrt((x  —  data$x[t])^2  +  (y  —  data$y [ t ] ) ^2) 

Pt  <-  1 

for(r  in  4: ncol ( data ) )  { 

rho.tsr  <—  sqrt  (  data  [  t  ,  r  ]  *  d.ts) 
switch  (  sensor  , 

cookie  =  if  (rho.tsr  <=  rho.O)  { Pt  <—  0;  break}, 

fermi  =  Pt  <—  Pt  *  (1  —  1  /  (1  +  1 0^(((  rho  .  tsr  /  rho.O)  —  1)  /  b))), 
exp  =  Pt  <-  Pt  (1  -  10^(-0. 30103  *  rho.tsr  /  rho.O)), 

stop(pasteO(’distribution  distr,  ’"  not  known  ’  ) ) 

) 

) 

return  (1  — Pt) 

) 
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lower .  bound  <—  function  ( data  ,  gamma,  rho_0  ,  obj  ,  sensor,  b,  gap_method)  { 

#  calculates  a  lower  hound  for  sector  gamma 

# 

i  Args  : 

#  data:  Data  set  created  hy  function  .  create .  data () . 

M  gamma:  Sector  data  created  in  function  DiBSO 

i  rho_0:  Range  of  the  day. 

M  long est _  edge :  Termination  condition, 

i  ohj:  Objective.  Possible  values:  *avg\  *min\ 

#  sensor:  Sensor  model.  Possible  values:  ^cookie\  ^fermi\  ^  exp  \ 
i  b:  Diffusivity  parameter  for  Fermi  model. 

#  gap_method:  Optimality  gap  method.  Possible  values:  ’none*,  ’center’,  ’corners’,  ’ 
corners  .  center  ’ . 

^  Returns : 

^  Lower  bound  for  sector  gamma. 
switch  ( gap_method  , 

none  =  return  (0)  , 

center  =  { 

potentials  <—  data. frame (x=mean( c (gamma$xmin  ,  gamma$xmax) )  , 

y=mean(  c  (gamma$ymin  ,  gamma$ymax) ) ) 


corners  =  { 

potentials  <—  data. frame (x=c  (gammaSxmin  ,  gammaSxmax  ,  gammaSxmin  ,  gammaSxmax)  , 

y=c  (gammaSymin  ,  gammaSymin  ,  gammaSymax  ,  gammaSymax)) 


center . corners  =  { 

potentials  <—  data,  frame  (x=c  (gammaSxmin  ,  gammaSxmax,  gammaSxmin,  gammaSxmax, 

mean( c (gammaSxmin  ,  gammaSxmax)))  , 
y=c  (gammaSymin ,  gammaSymin,  gammaSymax,  gammaSymax, 
mean(  c  ( gammaSymin  ,  gammaSymax) ) ) ) 

}, 

stop  (  pasteO  (’ gap  method  gap_method  ,  not  known’)) 

) 

potentials  Svalue  <—  supply  ( 1 :  nrow  (  p  o  t  e  n  t  i  al  s  )  ,  function(p)  { 
switch  (obj  ,  avg  =  value  <—  0,  min  =  value  <—  Inf, 

stop(pasteO(’objective  obj,  not  known  ’ ) ) ) 

if  (distr  ==  ’cookie’)  value  <—  0 
for(t  in  1 :  nrow(  data  ) )  { 

Pt  <—  .Pt(data,  t,  p  o  t  e  n  t  i  al  s  $x  [  p  ]  ,  p  o  t  e  nti  al  s  Sy  [  p  ]  ,  rho_0  ,  sensor,  b) 
if(sensor  ==  ’cookie’)  { 

value  <—  value  +  Pt*data[t,  ’value’] 

}  else  { 

switch  ( obj  , 

avg  =  value  <—  value  +  data[t,  ’value’]  *  Pt  /  nrow(data), 
min  =  value  <—  min(value  ,  data[t,  ’value’]  *  Pt) 

) 

} 
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return  (value) 


1) 

potentials  <—  p  o  te  nti  al  s  [  order(— p  o  te  n  t  i  a  1  s  $  value  )  ,  ] 
return  (  pot e  nti  al  s  $  value  [  1  ] ) 


#  Ex  amp  I  e 

#  The  next  examples  demonstrates  the  use  of  DiBS. 
targets  <—  data.franie( 


obj  ect 

=  c( 

’T1  ’  , 

’T2’  , 

’T3  ’  , 

’T4’  , 

’T5  ’  , 

’T6'  , 

’T7’  , 

’T8’  , 

’T9’  , 

’TIO’ ) 

X 

=  c( 

-15.00, 

6.25,- 

'  9.25, 

0.75,- 

-  4.00, 

9.00,- 

9.50  , 

2.75  , 

3.50,- 

4.00) 

y 

=  c(- 

-  4.75,- 

9.50,- 

-11.00, 

8.75  , 

5.50, 

-11.25,- 

15.00, 

11.50, 

-12.25  , 

7.50) 

value 

=  1) 

receivers 

< 

data,  frame  ( 

obj  ect 

=  c( 

’R1  ’  , 

’R2’  , 

’R3’) , 

X 

=  c( 

4.25  , 

14.25,- 

6.75), 

y 

=  c( 

-14.25, 

6.25  , 

12.50) ) 

#  1.  fermi  model,  maximize  average  detection  probability ,  longest  edge  <=  0.25 
sectors  <—  DiBS(targets  ,  receivers  ,  rho_0  =  3,  longest_edge  =  0.25) 

M  2.  cookie  cutter  model,  longest  edge  <=  0.1 

sectors  <—  DiBS(targets  ,  receivers  ,  rho_0  =  3,  longest_edge  =  0.1,  sensor= ’ cookie ’ ) 

#  3.  like  1.  but  longest  edge  <=  0.001,  use  center  for  lower  bound 

sectors  <—  sectors  <—  DiBS(targets  ,  receivers  ,  rho_0  =  3,  longest_edge  =  0.001, 

gap_method  =  ’center’,  gap  =  0.01) 


DiBS  Implementation  -  This  code  shows  a  basic  DiBS  implementation  in  R.  It  features  all 
sensor  models  mentioned  in  Section  1.2  as  well  as  longest  edge  and  optimality  gap  termination. 
It  assumes  that  plane  rotation  has  been  applied  prior  to  running  DiBS.  Uneven  rectangles  are 
handled  by  stopping  sector  splitting  in  the  direction  that  reaches  the  longest  edge  termination 
condition. 
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Supplementals 


•  Computer  Code  of  the  analysis  and  decision  tool  -  to  be  used  with  Appendix  A. 

•  Computer  Code  of  the  DiBS  algorithm  implemented  in  R  -  to  be  used  with  Appendix 

B. 

The  supplementals  are  available  at  Dudley  Knox  Library  of  the  Naval  Postgraduate  School 
in  Monterey,  CA. 
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